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ABSTRACT 


The  present  report  gives  a detailed  description  of  Program 
CABUOY.  which  analyzes  in  the  time  domain  the  two  dimen- 
sicnal  dynamic  behavior  of  genera!  ocean  cable  systems  consisting 
of  a surface  buoy,;  connecting  table,  and  intenrediatc  bodies. 

The  equations  which  model  the  motions  of  the  surface  waves 
and  the  various  components  of  the  cable  system  are  presented, 
.ind  the  subroutines  of  the  program  are  briefly  outlined, 
instructions  on  use  of  the  program  include  a listing  of  the  input 
KiiAD  .statements,  definitions  of  the  input  variables,  and  a 
number  of  comments  on  the  entering  of  input  data.  Several 
sample  problems  are  given  to  illustrate  use  of  the  program,  the 
output  ot  the  program,  and  computer  costs  lor  a range  of  cases. 
The  listing  of  the  ;;.i>gram  is  given  in  the  appendix.; 


AOMiNISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  authorized  by  the  Naval  Air  Development  Center 
under  Project  Orders  4-060!  and  0-0(ill  re.ipectivcly  dated  5 March  1974  and  1.^  February 
1976..  The  work  was  performed  under  interna!  Work  Units  M.S52-130  and  1-1552-145. 

INTRODUCTION 

The  dynamic  motion  characteristics  of  cabt;  system^  arc  cu'renlly  of  extreme  interest, 
and  several  major  surveys  of  cable  dynamics  studies  have  been  made  in  recent  years 

An  ocean  cable  system  generally  consists  of  the  following  three  components: 

1 . A ship  or  surface  float  at  the  upper  end . 

2.  A cable  whtise  properties  may  vary  along  its  length, 

3.  Intermediate  bodies  along  the  cable,  inc’.uding  the  possibility  of  a body  at  the  lower  end 
Previous  studies  have  usually  focused  on  only  one  of  the  above  components.  For  example, 
the  principal  emphasis  m many  studies  is  on  the  dynamic  characteristics  of  the  cable  ilsell. 


'DiOon.  D.B..  "An  Inventory  of  Current  Mathematical  Mtwleis  of  Scientiric  Data  - Gathering  Moors." 
Hydrospace-Challcnger.  Inc.  TR  4450  0001  (Feb  1973). 

A complete  list  of  references  is  given  on  pages  61-63. 

^Choo.  Y.I.  and  MJ.  Casarella.  "A  Survey  of  Aiulytical  Methods  for  Dynamic  Siniulation  of  Cable-Btidy 
Systems.”  Journal  of  iiydronaulics.  Vol.  7.  N«>.  4.  pp,  1 37  - 144  (Oct  1973). 

-^Aiberlsen,  N.D..  “A  .Survey  of  Techniques  for  the  Analysis  and  Design  of  Submerged  Moonng  Systems." 
Civil  Fngineermg  Laboratory  Technical  Report  R4tl5  I Aug  1974), 
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Conditions  at  the  ends  of  the  cable  are  then  either  those  of  prescribed  motions  or  simple 
representations  of  the  surface  buoy  or  lower  body.  It  appears  that  such  studies  were  earned 
out  mainly  to  demonstrate  the  feasibility  of  a panicular  method  of  solving  for  the  dynamic 
characteristics  of  the  cable..  Choo  and  Casarclla*  discuss  the  merits  and  drawbacks  of  the 
three  principal  analytical  methods:  linearized  frequency-domain  method,  method  of 
characteristics,  and  finite  element  method.  In  other  studies,  the  principal  emphasis  is  on  the 
dynamic  characteristics  of  the  surface  buoy  or  the  lower  body  and  the  effect  of  the  cable  is 
then  approximated  in  various  ways  It  is  clear  that  these  studies  are  suitable  only  for 
analyzing  particular  types  of  cable  systems,  also,  only  the  dynamic  characteristics  of  certain 
components  are  accurately  described. 

The  present  report  gives  details  on  Program  CABUOY.  which  analyzes  in  the  time  domain 
the  two-dimensional  dynamic  behavior  of  ail  three  components  of  a general  cable  system. 

This  program  has  already  been  briefly  described.^  Although  it  was  developed  principally  to 
analyze  the  dynamic  behavior  of  sonobuoy  systems,  for  which  it  is  of  interest  to  know  the 
dynamic  behavior  of  the  surfac'C  buoy,  connecting  cable,  and  lower  acoustic  detection  units, 
the  great  generality  a::d  versatility  of  the  program  make  it  useful  for  a wide  variety  of  other 
cable  systems. 

The  report  first  presents  in  detail  the  equations  which  form  the  basis  of  the  program. 
These  include  equations  for  the  steady-state  cable  configuration  and  for  the  dynamic'  motions 
of  the  surface  waves,  surface  buoy,  cable,  and  intermediate  bodies.  Eacit  of  the  program 
subroutines  is  briefly  descTibed.  Detailed  input  instructions  include  a listing  of  the  input 
FORTRAN  READ  statements,  definition  of  the  input  variables  contained  in  these  READ 
statements,  and  comments  on  the  entry  of  input  data.  Several  sample  problems  serve  to 
illustrate  program  use,  output,  and  computer  costs.  The  program  is  listed  in  the  appendix. 

Both  the  input  instructions  and  the  sample  problems  illustrate  the  wide  applicability  of 
the  program.  The  sample  problems  range  from  a parametric  study  of  the  accuracy  and 
computer  cost  of  various  finite  element  representations  of  the  cable  to  the  analysis  of  a 
complete  buoy-cable-body  system  moored  in  the  presence  of  typical  ocean  waves  and  current 
profiles.  The  characteristics  of  each  component  of  the  ocean  surfacx*  waves  may  be  sivcified 
by  the  user  or  may  be  internally  generated  by  the  program  by  means  of  the  Pu  rson-Moskowit/ 
spectrum.  The  surface  buoy  at  the  top  of  the  cable  may  have  a relatively  wide  range  of  si/es 
aitd  shapes.  It  may  be  a prolate  e>r  oblate  spheroid  of  any  aspect  ratio  provided  its  hon/ontal 

H.T..  “PreB.'ninary  Report  on  a Forlran  IV  Compu’er  Program  for  the  Two-DimensKinal  Dynamic 
Behavior  of  General  Ocean  Cable  Syiteim,'’  DTNSRDC  Departmental  Report  SPD-A.t.^^)!  (Aug  I *>75). 
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length  is  small  compared  to  the  wavelengths  of  the  significant  ocean  waves,  or  it  may  be  a 
spar  buoy  of  any  size..  The  reasons  for  these  particular  choices  of  buoy  shapes  and  sizes  are 
given  in  the  section  on  surface  buoys.  Alternatively,  motions  may  be  prescribed  at  the  upper 
end.  The  user  determines  the  accu;-acy  and  computer  cost  of  the  dynamic  analysis  of  the 
cable  by  specifying  the  total  number  of  cable  segments  as  well  as  the  length  of  each  segment. 
Several  different  formulations  are  given  for  the  added  masses  and  drag  coefficients  of  the 
intermediate  bodies. 


STEADY -STATE  CALCULATIONS 


CABLE  EQUATIONS 

The  program  first  calculates  the  configuration  of  the  cable  system  in  the  presence  of  a 
steady-state  cunent*  alone,  in  the  absenc'e  of  any  time-dependent  excitations.  The  differential 
equations  for  the  steady-state  configuration  of  the  cable  are  well  known  and  take  the  following 
form  for  the  coordinate  system  shown  in  Figure  1 ; for  example,  see  Springston.^ 


Figure  I — Definition  of  Coordinate  System 


*The  term  current  is  used  to  denote  the  steady-state  fluid  velocity  relative  to  the  cabla.  For  the  case  of  a 
cable  towed  in  still  water,  the  cunent  has  magnitude  equal  to  and  direction  opposite  to  the  towing  velocity. 

‘‘Springston,  G.B.  Jr..  “Generalized  Hydrodynamic  Loading  Functions  for  Bare  and  Faired  Cables  in  Two- 
Dimensional  Steady-State  Cable  Configurations.”  NSRDC  Report  7424  (Jun  1967). 
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■T^  + ! + sin0W  = O 


+ G + cos  0 W = 0 


(2) 


|i=,.e 


(3) 


dx  . 

( 1 + e ) sin  ^ 


(4) 


= ( 1 + e ) cos  d 

ds,. 


(5) 


where  T = cable  tension 

0 = angle  of  the  cable  segment  with  the  vertical  (see  Figure  I ) 

Sy  = reference  cable  length  (when  T = measured  from  the  initial  point 
= reference  tension 

I.G  = normal  and  tangential  drag  forces,  re.^pectbely.  per  unit  length  acting  on 
the  canle 

W = weight  in  fluid  per  unit  length  of  the  cable 
s = stretched  cable  length  measured  from  the  initial  point 
e = cable  strain;  e = 0 at  T = 

\ = horizontal  displacement,  positive  to  the  right 

y = vertical  displacement,  positive  downward 

For  smooth,  approximately  round  cables,  the  normal  and  tangential  drags  may  be  taken 
as  respectively  proportional  to  the  squares  of  the  velocities  normal  and  tangential  to  the  cable- 


I =-pCp  dCnk'n 


<(i) 


G = - p r,  d c,  I c, 


(7) 


where  p = fluid  density 

= normal  and  tangential  drag  coefficients  respectively 
d = cable  diameter 

c^  = compKineni  of  the  current  normal  to  the  cable  = c cos  O 

c,  = component  of  the  current  tangential  to  the  cable  = - c sin  0 

c = magnitude  of  the  current,  taken  to  act  on!\  in  the  \-direction 
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The  tension-strain  function  is  assumed  to  be  of  the  form^ 


where  C'|  = constant  of  elasticity:  Cj  - AF.  for  a linearly  elastic  cable 
A = cross-sectional  area  = ird*/4  for  a round  cable 
E = modulus  of  elasticity 

C,  = an  exponent;  C,  = 1 for  a linearly  elastic  cable 


(8) 


Equation  (8)  enables  a nonlinear  tension-strain  relation  to  be  modeled  by  only  two  input 
variables.  and  C,.  It  is  more  convenient  to  express  e as  a function  of  (T  -T,, ) in  order  to 
eliminate  it  in  Equations  (3»  through  i5>: 


€ = 


(9» 


INTERMEDIATE  BODIES 


It  is  assumed  that  conditions  at  the  top  of  the  cable  are  known  and  the  integration  of 
Equations  1 1 ' thtough  (5)  proceeds  down  the  cable.  .At  an  intermediate  body,  the  integration 
must  be  interrupted  and  the  unknown  cable  variables  T^^  and  0^  below  the  body  must  be 
related  to  the  known  vanabies  Ij.  and  0j.  above  the  body.  Reducing  <nc  three-dimensional 
et|uations  contained  in  References  7 ;,nd  8 to  the  present  twivdimensional  case  results  in  the 
following  two  e«iuations  for  T^  and 

Tu  = (sin  p|^  Tj^  \ B 


0 


I \ ^ 

cos  Pj.  T^  + Wg 


III) 


where  1)^  is  the  drag  on  the  Inidy  and 


^’Thresher.  R.W.  and  J.ll.  Nath.  ”Anchor-!jsS  Deployment  Procedure  for  Mooring."  Oregtm  Stale 
University  Report  7.^-5  (Jun  1973). 

^Wang.  H.T..  ■‘^.ffcct  of  Nonplanar  Current  Profiles  on  the  Configuration  of  Mtwred  Cable  Systems." 
NSRDC  Report  3692  (Oct  1971). 

^Wang.  M.T..  “A  FORTRAN  IV  Program  for  the  Tlircc-Dimcnsional  Steady-State  Configuration  of 
l-.xtcnsiblc  Flexible  Cable  Systems."  NSRFX'  Report  4384  (Sep  1974). 


.1:. 

Here  alM).  ^ |,  A^  is  ihe  drag  area  ot  the  body  tor  flow  in  tho  \-direction  and  ts  thi-  wi-ight 
of  the  bud.  'n  fluid. 


BOUNDARY  CONDITIONS 
Initial  Value  Cases 

The  integration  of  the  above  differential  equations  is  most  convenient  when  the  tension 
I and  the  angle  0 are  known  at  one  end  of  the  cable  I'liese  are  known  for  certain  cases  of 
single-point  moored  cables  and  towing  cables  for  the  moored  cases,  the  program  starts  with 
the  known  conditions  at  the  top  of  the  cable  and  integrates  the  five  differentia!  equatitms 
until  the  lower  end  ol  the  cable  is  reached.  This  js  the  simplest  case  for  the  program  since 
the  mimlvering  of  the  cable  segments  and  intermediate  bodies  starts  at  tl’.;  top  of  tlic  cable 
For  towing  cable  cases,  where  the  conditions  are  known  at  the  lower  towed  binlv.  the  program 
integrates  the  differential  equations  twice.  The>  are  first  integrated  from  the  towed  bodv  to 
the  upper  point,  thus  fixing  the  conditions  at  this  |H)int.  Then,  in  order  to  conform  to  the 
numbering  system  which  is  used  lor  the  dynamic  calculations,  the  equations  are  integrated 
on;e  again  from  the  upper  point  down  to  the  lower  lowed  body. 


BourKJary  Value  Cases 


In  main  applications  the  values  of  I and  0 are  not  known  a prion  at  any  point  along 
the  cable.  It  is  necessary  to  treat  these  cases  as  boundary  value  problems  .md  me  iteration 
techniques  to  oPlain  the  solution.  I he  present  program  contains  the  iteration  ^c!le^K's  l«ir 
two  c.ises  of  particular  interest  for  or.obiioy  systems  a cable  ol  given  length  moored  in  a 
given  ocean  depth''  .iiul  a tree-tlo.it ing  cable  system.'"  Ihe  iteration  subroutine  is  written  so 
th.it  the  User  m.iy  c<Mvciiientiy  imp’-, ment  iteration  schemes  for  other  applicMiioiw. 


'Uang.  M T and  B.I..  Webstci.  "Current  Prolilcv  Which  (iivc  Rise  to  Nonunniu:  Solutions  ol  Moored 
C'ahlc-  Systems.”  Paper  OTC  Fourih  Annual  Ottshorc  Technology  Conference.  Houston.  Tevas 

(May  i')7:). 

'"W'ang.  Il.T  and  T.I..  Moran.  ".Analysis  ol  the  Tw**-l)imcnsional  Sieady-Sialc  Behavior  of  I .xiensibic 
l tce-1'lo.il mg  Cable  Sysleins."  NSRIX'  Report  ,^7Z1  (Ocl  l'>71) 
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OCEAN  SURFACE  WAVES 


DESCRIPTION  OF  MOTION 


For  ocean  depths  greater  than  one-half  the  wavelength."  the  water  particle  trajectories 
due  to  a single  progressive  wave  are,  according  to  linearized  first-order  theor\',  given  by"  ’’ 


=a^ 


e“*‘y  cos  (kx  - ot  + * 


< I .^a  I 


- - e'^^y  sin  ( k\  - ot  + 0^. ) 


(l.sh) 


where  ~ particle  dispiaceirients  in  the  (x.yl  directions,  respectively 

a^  = wave  amplitude 

k - wave  number  = 2ir/X 

X = wavelength 

a = circular  frequency  = ^ 2»gA  = 2»f 
f = frequency 

t = time 

g = gravity  constant  = 32.2  ft/sec^  t‘).8l  m sec* » 

6^  = phase  angle 

It  is  of  interest  to  note  that  the  trajectories  describe  circular  orbits  with  a radius  which  decays 
exponentially  with  depth. 

For  an  irregular  sea  consisting  of  N distinct  components,  the  resultant  water  particle 
displacements  are  obtained  by  a summation  of  the  above  expressions,  resulting  in 


X 


w 


= Z)  3wi«'  COS  tkjX-0,t  +0^,1 
1 = I 


ll4.it 


k 

= sin(kjX-o,t +0^,t  «l4bi 

i = l 

Differentiations  with  respect  to  time  yield  the  following  results  for  water  particle  velocities 
and  accelerations 


’ * Lamb.  H.,  “Hydrodynamics.”  Sixth  Edition.  Dover  Publicatums.  New  York  1 1‘)45  j.  pp.  .th.t  - 170. 
pp.  152-  ;S5. 

*^Wehau$en.  J.V.  and  E.V.  Laitone.  "Surface  Waves.”  in  “Handbuch  der  Phssik.”  Vol  n.  Springer  Vcrlag. 
Berlin  (1^60).  pp.  446  -778. 
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(I  5a) 


= Z/  Oi  ^Wi  ' S'"  - °i*  + ^Wi  > 


-k  y 

y = o a - c ' cos  (k-x  - o t + 0 • ) 

■'W  1 WJ  ..w.i  Wji  wi 

1 = 1 


(15b) 


^ -ky 

'w  = Z ^ ' ‘■'OS  < kjX  - Ojt  ■^  e^j ) 

i = 1 


: I6a) 


^ > *k,y 

Vw  = Z Of  V '■•  S'"  - 0,‘  + ®wi‘ 

i = I 


(16b) 


CHOICE  OF  COMPONENTS 

The  computer  program  allows  the  user  two  options  for  describing  an  irregular  sea.  He 
may  sjKcify  the  values  of  N,  a^,.  Oj.  and  or  he  may  use  an  eneigy  spectrum  S^to)  to 
define  the  wave  amplitudes 

a^j  = y/  Sjlo,)  Ao  ( 1 7) 

where  Oj  - Oj  _ j + 

Ao  = (Oy  - o^)'N 

Oy  = upper  limit  of  the  significant  range  of  o’s 
Oy-  = lower  limit  of  the  significant  range  of  o's 

The  program  uses  the  Pierson-.Moskowitr  energy  sea  spectrum  of  the  form 

Sjto)  =— (18) 
o' 

where  A = O.OOS!  g-  and  B = 33. 5n  h,*,.  Here  h,,j  is  the  significant  wave  heign.t.  the  average 
of  the  one-third  highest  |H.*ak-ti>-t rough  heights.  .-\s  reported  by  Frank  .’-nd  Salvescn.''  the 
I Ith  International  Towing  Tank  Conference  (Tokyo)  ( 16(>p)  recommended  the  spectrum  in 
this  form  for  computations  when  information  is  not  available  on  typical  sea  spectra. 


rank.  W.  and  N.  5>alvcscn.  “The  frank  Oosc-l'ii  Ship-Motion  Computer  Prvigram.’*  NSRIXT  Reptut 
3:8*)  (Jun  l‘)70). 


S 


Other  forms  for  S^to)  can.  of  course,  be  conveniently  programmed.  Values  for  h,,j.  o^,  and 
Oj,  for  State  0 to  seas  may  be  found  in  Table  1 of  Reference  1.^.  The  program  sets  the 
values  of  0^1  to  be  evenly  spaced  from  to  .^60-d^,i  degrees. 

PRESCRIBED  SURFACE  MOTIONS 

The  program  allows  the  user  either  to  prescribe  the  motion  at  the  surface  or  to  describe 
It  by  means  of  differential  equations  ot  motion  for  i surface  buoy 

If  the  surface  buoy  or  ship  is  sufficiently  large  that  its  motions  are  not  appreciably 
affected  by  the  presence  of  the  cable,  these  motions  may  be  calculated  separately  and  used  as 
input  for  the  present  program.  Several  programs  are  available  to  calculate  the  pitch  and  heave 
motion  responses  of  surface  ships,  for  example,  the  Frank  (lose-Fit  Ship  Motion  Computer 
Program.'^  This  approach  is  also  valid  for  laboratory  simulations  of  cable  dynamics  where 
the  motions  at  the  upper  end  of  the  cable  are  often  prescribed.  third  area  of  application 
could  be  a full-scale  trial  where  the  motions  of  the  surface  ship  or  platform  can  he  readily 
measured. 

The  program  considers  the  prescribed  motion  of  the  end  of  the  cable  as  composed  of  a 
series  of  sinusoidal  components  in  the  horizontal  and  vertical  directions 

N 

Xj  = ^ a^j  cos  (-  :jrf|t  + 1 1‘>) 

I = I 


N 

Vj  = -Uy,  sin(-:irf,t (20) 
i = I 

where  x^.  y^  = horizontal  and  veitical  components  of  the  surface  motion,  respectively 
^xr‘*vi  ~ amplitudes  for  the  ith  component  of  and  y^.  resps*cti\ely 
f|.  Oji  = frequenev  and  phase  angle  for  the  ith  component,  respectivelv 

If  so  desired,  other  forms  for  the  prescribed  motion  may  be  comeniently  added  to  the 
program,  e.g..  other  functions  of  time  such  as  powers  ol  t or  exponentials 


SURFACE  BUOY  EQUATIONS 


GENERAL  CONSIDERATIONS 

it  is  well  known  that  the  added  muss  and  damping  coefficients  of  surface  buoys  are.  in 
general,  functions  of  the  frequency  of  the  oscilbtion.*^  In  the  time  domain,  this  requires 
the  solution  of  integrodiffercntial  equations  which  contain  convolution  integrals.  Alternatively, 
if  the  frequency-dependent  coefficients  can  be  expressed  as  simple  polynomials  of  the  freque.ncy, 
the  integrodiffercntial  equations  may  be  replaced  by  a set  of  higher  order  differential 
equations.*'*  In  either  case,  the  solutions  are  complex  and/or  time-consuming  in  the  time 
domain,  llius,  surface  buoy  motions  have  usually  been  solved  in  the  frequency  domain,  in 
this  approach,  the  steady-state  harmonic  response  is  obtained  for  each  frequency  component 
of  the  exciting  surface  waves.  The  total  response  to  the  sum  of  the  individual  wave  compo- 
nents is  then  obtained  by  linear  superposition.  Experiments  have  shown  that  this  procedure 
generally  yields  satisfactory  results  for  pitch  and  heave  motions  of  surface  ships. 

Because  of  difficulties  in  solving  buoy  equations  in  the  time  domain,  the  frequency 
domain  approach  has  also  been  used  to  study  the  motion  of  cable-buoy  systems.  Perhaps  the 
most  comprehensive  of  these  studies  is  the  Goodman  et  al.  computer  program'^  which 
considers  four  different  buoy  shapes.  In  addition  to  facilitating  the  solution  for  general  buoy 
shapes,  the  frequency  domain  approach  has  the  additional  advantages  of  immediately  giving 
the  steady-state  harmonic  response  (no  need  to  wait  for  the  transient  response  to  die  down) 
and  of  reducing  the  computer  time  required  to  obtain  cable  motions.*’  However,  the  draw- 
backs to  this  approach  include  neglect  of  all  nonlinearities  and  the  assumption  that  ail  the 
dynamic  response  variables  are  small  compared  to  their  steady-state  values.  This  approach 
would  not  be  able  to  predict,  for  example,  the  large  dynamic  snap  loads  which  occur  when 
the  cable  goes  slack. 

In  view  of  the  above  drawbacks  and  also  view  of  the  existence  of  the  comprehensive 
freo.uency-domain  computer  program  described  in  Goodman  et  al..*^  it  was  decided  to  use  a 


’■*Tick.  I..J..  “Difcrcniul  Equations  with  Frequency-Dependent  Coefficients.”  Journal  of  Ship  Research. 
Vol.  3.  No.  2.  pp.  45  -40  (Oct  1959). 

"Ogilvie.  T.F..  “Recent  Progress  toward  the  Understanding  and  Prediction  of  Ship  Motions.”  Fifth 
Symposium  of  Naval  llydrodyrumics.  Bergen.  Norway,  pp.  3-128  (Sep  1964). 

**’(HH>dman.  T.R.  et  al..  "Static  and  Dynamic  Analysis  of  a Moored  Buoy  System.”  National  Data  Buoy 
Center  Report  6113.1  {.4pr  1972). 

*’Wang.  H.T..  “A  Two  Dcgree-of-Frcedom  Mode)  for  the  Two-Dimensional  Dynamic  Motions  of  Suspended 
Extensible  Cable  Systems.”  NSRIX"  Report  .3663  (Oct  1971). 
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lime  domain  approach  in  the  present  study.  In  order  to  make  this  approach  leasible.  it  was 
important  to  find  classes  of  buoys  which  did  not  have  frequency-dependent  added  mass  and 
damping  coefficients.  A literature  search  revealed  two  such  classes'  spar  buoys®*-*’  and  small 
buoys.  Spar  buoys  are  buoys  with  circular  cross  sections  and  large  draft-to-diametcr  ratios 
Hi'b.  Because  of  their  slenderness,  the  added  inertia  terms  of  these  buoys  are  essentially  those 
for  infinite  fluid,  and  the  frequency-dependent  wave  damping  coefficients  are  zero  to  first 
order  approximation.**  The  other  class  corresponds  to  buoys  whose  typical  dimension  a is  so 
small  compared  to  the  ocean  wavelengths  X that  the  reduced  frequency  a given  by 

o = 2rr^«l  tell 

is  much  los;  than  unity  for  the  range  of  X values  corresponding  to  ocean  waves  of  m'.erc'-^t. 

For  surfac*e  buoys  of  sonobuoy  systems,  whose  typical  dimension  is  ol  the  order  of  1 ft 
(0.305  m).  the  above  condition  holds  for  the  large  majority  of  sea  states.  When  hquation  (21 ) 
holds,  the  wave  damping  terms  go  to  zero  and  the  added  inertia  terms  for  o = 0 may  be  used. 
In  this  case,  the  ocean  surface  behaves  essentially  as  a rigid  plane.*"  and  tor  the  case  of  a buoy 
whc>se  axis  of  symmeiri'  is  aligned  with  the  y-axis.  the  added  mass  m surge  is  equal  to  the 
infinite  fluid  value.  The  added  mass  coefficients  for  pitch  and  heave  must  l>e  calculated 
separately  for  each  shape  considered.  Since  thes**  cwfllcients  have  Ix’en  studied  for  several 
cases  of  jblate  and  prolate  spheroids.*' and  also  because  they  represent  mathematical 
shapes  which  are  similar  to  surface  buoys  of  sonobuoy  systems,  it  was  decided  to  represent 
the  small  buoys  by  prolate  and  oblate  spheroids.  Both  types  of  spheroids  .ire  characterized  b\ 
having  two  of  their  throe  axes  equal  in  length.  The  limiting  cases  for  a prolate  spheroid  are  a 


**Ncwm3n.  J2x'..  “The  MoIkmu  of  a Spar  Buoy  in  Regular  Waves.”  David  Taylor  M-.>dcl  Basin  Report 
1409  (May  1963). 

* 'Rudnick.  P...  “Motion  of  A Large  Spar  Buoy  in  Sea  Waves.”  Jourtal  ol  Ship  Research.  Vol.  1 1.  No.  4. 
pp.  25‘7. 267  (Dec  1967). 

•"Newman.  J.N...  "Marine  Hydrodynamics  (Ixcsure  Notes).’’  M.l.T  Depi  Na>.  .Arch,  and  Mar  i ng. 

(Spring  Term  1971 ). 

* ' Havelock.  T..  “Waves  due  to  a Floating  Sphere  Making  Periodic  Heasm^  Osciilaiions.”  Proceedings  ol 
!he  Koval  Society.  Vol.  231.  Serws  A.  pp.  1-7  (Ju!  J9.SSj. 

■*MacCamy.  R.C.>  “On  the  Heaving  Motion  of  Cylinders  ol  Shallow  Drali.”  Joninai  ol  Ship  Researcli. 
V'ol.  ' No  ,3.  pp  .34-43  iDcc  l*^6t). 

‘ ■'Kim.  W’.D  , “On  ihc  Forced  Oscillitions  of  Shaliow-Drali  Ships.”  Journal  of  Shir  Research.  Vv'l.  7. 

No.  2.  pp.  7-i.S  (Oc!  t‘»63). 

'■*Kim.  W'  I) . "f)n  s r»cc-noating  Ship  in  Waves.”  .iivurnal  tii  Ship  Research.  Voi  I()  No  .3.  pp  !S2  • *‘*1 
(Sop  I'3n>) 
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long  thin  cylinder  and  a sphere,  and  the  limiting  cases  for  an  oblate  spheroii  are  a thin 
circular  disk  and  again  a sphere.  It  can  be  seen  from  these  limiting  cases,  that  prolate  and 
oblate  spheroids  can  be  used  to  generate  a wide  range  of  shapes. 

DIFFERENTIAL  EQUATIONS  OF  MOTION 

The  linearized  differential  equations  for  the  surge,  heave,  and  pitch  motions  of  a surface 
buoy  or  ship  freely  floating  in  an  inviscid  fluid  are  well  known;  see.  for  example,  Frank  and 
•Salvesen'^  and  Newman.'*  These  equations  usually  contain  the  inertia  forces  (including  the 
added  hydrodynamic  inertia  forcesl.  the  wave-damping  forces,  the  exciting  forces  due  to  the 
incoming  ocean  waves,  and  the  restoring  forces  due  to  buoyancy.  The  lesulting  equations  are 
usually  solved  in  the  frequency  domain.  As  mentioned  previously,  the  wave-damping  forces 
may  lx*  neglected,  to  first  order,  for  the  two  types  of  buoys  considered  in  the  present  report., 
i he  fu'ces  due  to  viscous  drag  and  cable  tension,  which  are  not  usually  considered  in  the 
.dr.n-:  studies,  are  included  in  the  present  formulation.  The  viscous  drag  forces  which  are 
quadratic  in  the  motion  velocities  of  the  float,  are  usually  omitted  since  they  make  the 
equations  nonlinc.ir  and  complicate  solutions  in  the  frequency  domain.  The  inclusion  of 
these  forces  poses  no  problem  in  a time-domain  analysis.  The  cable  forces  arc.  of  course, 
zero  for  a freely  floating  buoy. 

If  the  above  forces  acting  on  the  buoy  (shown  in  Figure  21  are  considered  and  the  pitch 
angle  is  taken  to  be  small  such  that 

sin  i|/  =55  ^ ( 22a( 


cos  =^  1.  (22b) 

the  threi  differential  equa-.ions  for  the  surge  heave  f.  and  pitch  \!/  are  as  follows 

= (m  + P vi]  i + 1 A^^  = - f f (y  - y,.  ) kjty)  Sty ) dyl  it  = FK^  + l\  + + T,^^ 

^ (2.1) 

[%,  = ""  * S„pV|]  f = -l>>.  -T„  C4, 


0 = pgVBGC/  + FK. 


' ry  + r, ) Ty  - (r^  C'  + Ty ) - r^y  T, 


12 


where  A^.  and  A^^  = inertia  coefficients  defined  above 

m = mass  of  the  buoy 

Kg  and  Kjj  = added  mass  coefficients  for  surge  and  heave,  respectively 
V = submerged  volume 

H = draft 

y^.  = distance  of  the  center  of  gravity  below  the  undnturbed  free  surface 

kg(y»  = local  added  mass  c-.>efficient  for  surge 

Sty  I = local  cross-sectional  area 

FK^.  FKy  = exciting  forces  due  to  surface  waves  in  the  x-  and  y-directions. 
respectively 

= exciting  moment  about  the  center  of  gravity  due  to  surtace  waves 

= viscous  drag  forces  in  the  x-  and  y-directions.  respectively 

= components  of  the  cable  tension  in  the  x-  and  y-directions. 
respectively,  at  the  alt  ichment  point  to  the  buoy 

= wind  loading  on  the  buoy  in  the  x-direction 

= steady-state  component  of  tension  in  the  v-direction  at  the 
attachment  point  to  the  buoy 

= waterplane  area  of  the  buoy 

= vertical  displacement  of  the  ocean  surface,  defined  in  Equation  1 14b) 

= moment  of  inertia  anout  the  center  of  gravity 

= yB->'o 

= distance  of  the  center  of  buoyancy  below  the  undisturbed  free  surface 

= horizontal  and  vertical  distances  measured  from  the  center  of  gr.ivity 
to  the  cable  attachment  point 

■-  vertical  distance  measured  from  the  center  of  gravity  to  the  center  of 
the  v/ind  loading  force 

In  Equations  (2.^)-(25).  a dot  denotes  differentiation  with  respect  to  time. 

The  added  inertia  coefficiems.  presented  below,  are  all  calculated  for  the  case  v = 0 It 
should  be  noted  that  the  moment  terms  on  the  right-hand  side  of  Equation  i2.s)  may  K- 
readily  changed  to  account  for  large  values  of  which  negate  approximations  (22ai  and  i22hi. 
Mov.-ever.  the  added  inertia  terms  on  the  left-hand  side  of  Equations  l2Jt)  through  (2.^)  must 
K*  calculated  (by  pote  itial  flow  .ncthods)  for  each  new  value  of  C'.  For  example,  tor  .i  large- 
value  of  \1/.  the  added  mass  in  surge  for  a small  spheroidal  buoy  is  no  longer  equal  to  the 
infinite  fluid  value  since  in  this  case,  there  will  K'  motion  of  the  Iluid  perpendicular  to  the 
free  surface. 

Under  static  conditions,  the  submerge'^  volume  V must  support  both  the  weight  of  tiie 
buoy  in  air  (mg)  and  the  vertical  component  of  the  steady-state  tension  < l^^i 


FK^ 

Dx-Dy 

T.S  Ty 


T 

*twx 

T 
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yw 
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(26a) 


pgV  = mg  + Tyj 


P Pg 


(26b) 


If  the  input  dimensions  for  the  draft  and  cross-sectional  areas  of  tlie  buoy  are  such  that  the 
submerged  volume  does  not  equal  the  value  given  in  Equation  (26b).  the  program  internally 
multiplies  the  cross-scctional  areas  by  a constant  factor  so  that  the  submerged  volume  becomes 
exactly  equal  to  the  value  given  by  thi:.  equation.  .Ml  of  the  added  inertia  and  wave-e.xciting 
forces  given  in  the  following  sections  are  based  on  this  volume. 

The  lollowing  two  sections  present  derivations  for  the  added  inertia  coefficients  K^.  Kjj. 
and  the  wave-exciting  forces  FK^,,  FK^  for  the  two  classes  of  buoys  considered: 
spar  buoys  and  small  spheroidal  buoys. 


SPAR  BUOYS 

Because  of  the  slenderness  of  the  spar  buoy.  Newman'**  shows  that  its  added  inertias  for 
surge  and  pitch  are  identical  to  those  in  infinite  fluid.  Since  the  spar  buoy  has  a circular 
cross  section,  the  surge  inertia  coefficients  kj  and  K^.  arc  both  equal  to  I.  leading  to  the 
following  definitions  for  the  added  surge  and  pitch  inertia  terms  in  Equations  (23)  through 
(2.*:.: 

K3=k5=l  (27a) 

m + K3pV  = ni  + pV  (27b) 


>/ 


H 

(y->j.  )S(y)dy 


(27c) 


p f (y k5(y)S(y)dy  = 


M 

<y-Vf;)-  S(y)dy 


(27d) 


Newman  takes  the  buoy  to  be  sufficiently  slender  so  that  the  added  mass  for  heave  may 
be  neglected.  Adec  and  Bai**  have  shown  experimentally  that  for  the  case  of  a circular 
cylinder,  it  is  more  accurate  to  add  a term  corresf.jnding  to  one-half  the  added  mass  of  a 
circular  disk  (with  the  same  diameter  as  that  of  the  cylinder)  heaving  in  infinite  fiuid. 


'•''.Adce.  B.H.  and  KJ.  Bai.  ‘'Experimental  Studies  of  the  Behavior  of  Spat  Type  Stable  Platfcrms  in  Waves.” 
University  of  California  (Berkeley)  Report  NA-70-4  (Jul  1^‘70). 
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This  correction  term  has  been  incorporated  into  the  present  formulation  by  takir  > the  radius 
of  the  disk  to  be  the  mean  radius  of  the  spar  buoy  r.  defined  by 


jrr*  H = V.  f = 


Since  one-half  of  the  added  inertia  of  a circular  disk  heaving  m infinite  fluid*'  is  (4/3)  pr^, 
the  following  expression  is  obtained  for  Kjj  in  Equation  (24) 


K„ 


i i-3 

3 V ^ ^ 3 V [uHl 


(28) 


Newman  shows  that  the  wave-exciting  forces  are  simply  the  Froude-Krylcv  forces,  which 
may  be  obtained  by  integrating  the  pressure  field  generated  by  the  ocean  waves  around  the 
contour  of  the  buoy.  The  following  expressions  arc  obtained  for  FK,^  and  FK^ 


where 


FK 


.X 


N 

- - Z COS  (kjX  - Oji  + Oo(kj)] 

i = i ^ 

N 

- Z [0|*  »wi  COS  (k,x  - Ojt  + 0,  (ki] 

i = l ^ 


Z'**  -k.y 

Qo^kj)  = p / e ' S(y)  dy 
''0 


Q,(kj)  = p 


(y-y,. ) Sty)  dy 


(29) 


(30) 


For  the  heave  motion,  the  Froude-Krylov  force  given  by  Newman  has  been  modified  to 
account  for  the  additional  heave  added-mass  term  given  in  Equation  (28).  This  modification 
has  been  so  made  that  in  the  limiting  case  of  a very  small  spar  buoy  (which  would  follow  the 
motions  of  the  waves  in  the  absence  of  cable  forces),  the  wave-exciting  term  in  Equation  (24) 
would  be  exactly  equal  to  the  inertia  term  on  the  left-hand  side  of  the  equation.  The  resulting 
force  has  the  from 


FKy  = [l  + I P **w,  ‘‘i"  < k,x  - Ojt  + 0^, ) 0„(  k, )] 


(31) 


The  expression  given  by  Newman  does  not  contain  the  correction  term  (4/3)p(V/jrH)^'*. 
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SMALL  SPHEROIDAL  BUOYS 


As  mentioned  previously,  for  the  limiting  case  of  ?cro  lediiced  frequency  where  the  free 
surface  behaves  as  a rigid  plane,  the  coefficient  for  the  surge  added  mass  Kj  is  equal  to  the 
infinite  fluid  value.  These  coefficients  may  be  calculated  by  the  formulas  and  tables  given  in 
Lamb." 

The  added  inertia  coefficients  for  heave,  pitch,  and  coupled  pitch-surge  motions  which 
have  components  normal  to  the  free  surface  are  not  the  same  as  the  coefficients  for  infinite 
fluid  Instead,  they  must  be  calculated  separately  to  incorporate  the  rigid  free-surface  ^ 
condition.  Inspection  of  Equations  (23)  and  (25)  shows  that  the  added  inertia  terms  for  ' 
pitch  and  coupled  pitch-surge  motions  may  be  defined  as  follows 


H 

(y  - ks<y>  sty)  dy  = 


■*■>0  * k5(y)S(y)dy 


= V - 2y^  (-  ) r^.  V r y^.'  Kj  \ (32a) 

M 

‘y  ~ 

where  r^  is  the  maximum  radius  of  the  buoy  and  is  used  to  render  the  coefficients  pjj  and 
>1,5  dimensionless. 

Thus,  a calculation  of  Pjj.  >r,j.  and  K„ . along  with  the  values  of  Kj,  for  infinite  fluid  as 
given  in  Lamb,  completely  determines  all  the  added  inertia  terms  in  Equations  (23)  through 
(25).  Bai*^  has  used  a finite  element  approach  to  calculate  coefficients  /ijj.  >i,j.  and  K„  for 
spheroids  with  draft  to  maximum  radius  ratios  (H/r^. ) ranging  from  0.1  to  10.  His  results 
agree  well  with  previous  results*'  *^  for  corresponding  cases  at  the  zero  reduced  frequency 
limit. 

Principally  for  the  sake  of  programming  ease,  the  present  program  consio^’s  only  the 
results  for  the  case  where  the  maximum  radius  lies  at  the  free  surface.  Consideration  of  other 
radii  at  the  free  surface  would  introduce  a-lditional  parameters  to  a description  of  the 
submerged  buoy.  In  addition,  it  is  expected  that  in  most  cases  the  maximum  radius  will  be 
close  to  the  free  surface  because  most  of  the  volume  (and  hence  buoyancy  capability)  of  the 
spheroid  is  concentrated  in  the  region  around  the  maximum  radius.  Results  for  several  cases 
where  the  waterline  does  not  occur  at  the  maximum  radius  arc  presented  in  Bai.’* 

^^Bai.  KJ..  “The  Zero-Frequency  Hydrodynamic  Coeffiaents  of  Vertical  Axisymmetric  Bi>difS  at  a Free 
Surface.”  Journal  of  Hydronauticc  (Jan  1^77). 
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In  the  range  0.1  < H/r^  < 10.,  the  coefficients  ^15.  #^55,  Kjj,  and  IC5  are  obtained  by 
linearly  interpolating  between  the  values  shown  in  Table  1 . 


TABLE  1 - VALUES  OF  ADDED  INERTIA  COEFFICIENTS 
FOR  SPHEROIDAL  BUOYS  AT  ZERO  REDUCED  FREQUENCY 


H/r 

w 

1h 

0.1 

0.074 

12.84 

1.27 

0.3 

0.2 

0.143 

5.84 

0.5S 

0.264 

0.3 

0.20.. 

3.672 

0.312 

0.237 

0.5 

0.31 

2.005 

0.117 

0.177 

0.7 

0.397 

1.323 

0.0358 

0.11 

0.9 

0.469 

0.96 

0.0038 

0.0039 

1.0 

0.50 

0.836 

0. 

0. 

1.5 

0.622 

0.484 

0.0731 

-0  191 

2. 

0.704 

0.323 

0.272 

-0.391 

3. 

0.804 

0.18 

0.993 

-0.788 

5. 

0.894 

0.082 

3.65 

-1.565 

7. 

0.933 

0.049 

8.20 

-2.38 

10. 

0.96 

0.028 

20.0 

-3.75 

The  values  of  Kj  are  those  given  by  Lamb,**  and  the  values  of  /i,j.  Pjj,  and  are  the 
values  given  by  Bai.^  subject  to  two  modifications  at  H/r^  = 10.0.  The  values  of  #i|5  and 
/ijj  for  H/r^  = 10,  respectively  20.0  and  -3.75,  correspond  to  those  obtained  by  simply  using 
strip  theory  with  kjty)  = 1 for  0 < y < H.  The  corresponding  values  of  and  Pjj  obtained 
by  Bai  are  respectively  17.84  and  -3.53.  These  modifications  were  made  principally  for  the 
sake  of  providing  continuity  with  the  approximation  used  in  the  range  H/r^  > 10.0,  where 
the  buoy  is  treated  essentially  as  a spar  buoy.  There  are  two  reasons  for  the  differences 
between  the  strip  theory  and  theoretical  finite  element  calculations.^  First,  the  strip  theory 
approach  neglects  the  flow  around  the  lower  end  of  the  buoy,  where  k^ly)  < 1.  Second,  the 
finite  element  representation,  where  only  the  nodes  of  the  elements  are  on  the  surface  of  the 
buoy.,  effectively  models  a smaller  buoy.  Both  of  these  effects  serve  to  make  the  strip  theory 
values  higher  than  the  corresponding  finite  element  results. 

W«v«  Exciting  Foron 

The  exciting  forces  FK^^  and  FKy  in  Equations  (23)  and  (24)  were  written  in  a form 
which  assumes  that  the  buoy  follows  the  wave  motion  in  the  absence  of  cable  forces  and 
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coupling  between  pilch  and  surge  motions.  This  is  the  case  for  buoys  with  dimensions  which 
are  small  compared  to  the  lengths  of  the  exciting  waves.  The  resulting  equations  take  the  form 

N 

FK^  = (1  + KjlpVx^  = -(1  + KjlpV  [oj*  a^jCosi.k.x  - Ojt  (33) 

1-1 


N 

FKy  = ( 1 + K„  ) p Vy^  = ( 1 + K„  i pV  [o?  a^,  sin  ik,x  - o,t  + <34l 

I = ! 

V y 

The  term  e ' which  appears  in  Equations  (16a)  and  i I6bl  for  and  y^.  ha%  been  omitted 
in  the  above  equations  since  it  is  ^ I under  the  assumption  of  small  reduced  frequency. 
Equation  (21 ). 

The  pitch-exciting  moment  is  computed  by  noting  that  because  of  the  sjinmetry  of  the 
buoy  about  the  vertical  axis,  only  the  horizontal  wave  motions  make  a contribution  to  pitch, 
resulting  in 


H 

(y-V(;»\v  + ks(y)l  S(y)  dy 


= p(-BG  V + p,jr^V  + y^.  KsV)x^ 


= -pV  (-  BG  + p,5  V + y^  Ks  > £ [o,-  a^,  cos  (k,x  - o,t  + (35) 

I = 1 

-k  \ 

Again,  under  the  assumption  of  small  reduced  frequency,  the  term  e ''  has  been  set  equal  to 
1 in  the  expression  for  x^,. 


DYNAMIC  CABLE  EQUATIONS 
GENERAL  CONSIDERATIONS 

As  mentioned  previously,  equations  for  the  present  study  are  solved  in  the  time  domain. 
Previous  cable  studies  have  considered  two  major  approaches  m the  time  domain*  the  method 
of  characteristics  and  the  finite  element  method.  The  method  of  characteristics  is  an  elegant 
method  which  reduces  the  originai  set  of  partial  differential  equations  to  a set  of  ordinary 
differential  equations  which  are  integrated  along  characteristic  lines  or  wavefronts.  This 
method  furnishes  valuable  lasight  into  the  various  modes  cf  cable  motion  but  solution  times 
are  typically  very  large.* 


The  finite  clement  method  seeks  to  represent  the  actual  cable  system  by  a series  of 
segments  and  nodes.  The  original  set  of  partial  differential  equations  is  then  reduced  to  a set 
of  ordinary  differential  equations  of  motion  for  the  nodes.  This  method  facilitates  the 
modeling  of  nonuniform  properties  along  the  cable  as  well  as  the  presence  of  intermediate 
bodies.  The  method  is  also  quite  flexible  in  that  the  number  and  location  of  nodes  is  left  to 
the  judgment  of  the  user.  In  arriving  at  his  selection  of  nodes,  he  may  consider  such  factors 
as  the  type  and  accuracy  of  the  dynamic  information  desired,  the  amount  of  computer  time 
available,  the  complexity  of  the  cable  system,  and  the  spatial  variation  of  the  environmental 
velocity  profiles. 

Principally  for  reasons  of  generality  and  flexibility,  the  finite  element  approach  was  used 
io  mode!  the  cable.  Straight  elements  are  used  in  the  present  formulation.  Webster^’  studied 
the  use  of  higher  order  curved  elements  to  model  cable  shape  and  concluded  that  the  first 
order  straight  element  appears  to  be  the  most  cost-effective.  In  particular,  he  showed  that 
one  second-order  quadratic  element  would  have  to  be  as  accurate  as  at  least  eight  straight 
elements  before  its  use  would  be  economical. 

Before  deciding  on  the  final  formulation,  a number  of  preliminary  approaches  for 
obtaining  the  diflerential  equations  of  motion  for  the  nodes  were  explored.  Two  approaches 
In  particular  were  considea'd  in  some  detail. 


PRELIMINARY  APPROACHES 

In  one  approach,  the  equations  were  formulated  in  a coordinate  system  aligned  with  the 
cable  segment.  This  is  the  approach  used  by  Rupe  and  Thresher®  to  obtain  the  dynamic 
motions  of  an  incxtensible,  uniform  cable.  For  the  present  c~ase  of  an  extensible  cable,  the 
two  unknowns  are  the  inclination  and  stretch  of  each  segment.  This  is  the  most  natural  way 
of  describing  the  configuration  of  a segment.  In  addition,  certain  cable  forces  such  as  tension, 
added  inertia,  and  drag  forces  are  most  conveniently  expressed  in  directions  normal  and 
tangential  to  a cable  segment.  However,  in  the  presence  of  intermediate  bodies  along  the 
cable,  for  which  the  inertia  and  drag  arc  most  conveniently  expres,sed  in  the  spatial  x-  and 
y-directions.  the  resulting  equations  are  greatly  complicated  by  the  transformation  required  to 
exprevs  the  body  forces  in  the  cable  coordinate  system.  The  cable  system  considered  by  Rupe 
and  Thresher  is  free  of  intermediate  bodies. 

•^Webster.  R.L..  “An  ApphcalKin  of  the  Finite  Clement  Method  to  the  f ^termination  of  Nonlinear  Static 
and  Dynamic  Resprmscs  of  Underwater  Cable  Structures,**  (General  Electric  Report  R76F.MH2  (Jan  l‘J76>. 

* Rupe.  R.C  and  Thresher.  R.W.,  “The  Anchor-l,ast  Deployment  Problem  for  Inextensible  Mooring 
Lnes.”  /\SME  Paper  74-WA;’(K.T-5  (Dec  Id74l. 
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Wang^  has  shown  that  relatively  fe<«b.5(ginents  are  required  to  accurately  describe  the 
overall  steady*state  conflguration  of  a cable;  these  results  suggested  a second  novel  approach. 
The  cable  was  conventionally  divided  into  a number  of  straight  segments  and  two  differenli^ 


equations  in  the  x-  and  y-directions  were  written  for  the  nodes  at  the  ends  of  the  segments. 
However,  each  straight  segment  was  subdivided  into  a number  of  intermediate  nodes,  as 
shown  in  Figure  3. 


Since  tiiese  intermediate  nodes  were  forced  to  move  along  the  straight  cable  segment,  only 
one  differential  equation  was  needed  to  describe  their  longitudinal  motion.  The  principal 
intention  of  this  approach  was  to  have  the  end  nodes  describe  the  overall  cable  configuration 
and  the  intermediate  nodes  describe  the  variation  of  tension  along  a cable  segment.  There 
was  not  sufficient  time  in  the  present  study  to  fully  explore  this  approach.  However,  it  was 
found  that  a certain  amount  of  bookkeeping  was  required  in  the  program  to  differentiate 
between  the  “end”  and  "intermediate”  nodes.  Also,  although  this  approach  reduces  the  total 
number  of  differential  equations  from  that  required  by  more  conventional  approaches,  the 


H.T.,  “Determination  of  the  Accuracy  of  Segmented  Represenutiotu  of  Cable  Shape.”  Journal  of 
Engineeriiy  for  Industry.  Vol.  97,  No.  2,  Series  B.  pp.  472-478  (May  I97S). 
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D,  = tangential  drag  force  acting  on  the  cable 
Djj  = drag  force  acting  or  the  intermediate  body 
Wj,  = weight  of  the  cable  in  fluid 
Wg  = weight  of  the  intermediate  body  in  tluid 
M = total  number  of  nodes 


Detailed  definitions  of  the  cable  forces  are  now  given,  and  the  drag  and  inertia  forces  for  the 
intermediate  bodies  arc  defined  in  the  next  major  section. 

It  IS  most  convenient  to  solve  the  set  of  2M  differential  equations  i3b»  and  (.^7»  with  Xj 
and  y,  uncoupled,  as  follows 


Xj  = <J  F - Kj  F )/tJ  J - K:‘  I 


1 "1  yi 


(3«a) 

(3ob» 


DEFINITION  OF  CABLE  FORCES 


Inertia  Forces 


Each  cable  segment  is  taken  to  be  a long  thin  cylinder  for  which  fluid  inertia  is  added 
only  for  acceleration  normal  to  the  segment.  Thus,  in  a coordinate  system  aligned  with  the 
segment,  the  inertia  force  F,^-  for  the  cable  segment  is  simply  given  by 


F„.  =(pk„+apAk„»..„+|ii„a, 


I.VM 


where  n = mass  per  unit  length  of  the  cable 

= reference  length  of  the  cable  segment 
a = added  mass  coefficient  = 1 .0  for  a round  cable 
\ = cros,s-scctional  area  of  the  cable 

ajj.aj  = accelerations  respectively  normal  and  tangential  to  the  cable  segment 

In  this  equation,  an  arrow  denotes  a vector. 

In  the  fixed  x-  and  y-directions,  for  which  the  ditferential  equations  are  wiitten.  the 
inertia  coefficients  are  not  constants  but  rather  are  lunctions  of  cable  inclination  <i>.  For  the 
coordinate  systems  shown  in  Figure  4.  a^^  and  aj  are  related  to  \ and  y.  respectively  the 
accelerations  in  the  x-  and  y-directions,  by 


Ujj  = x’  cos  0 - y sin  0 
a,  = - X sin  0 + y cos  <t> 


t40a» 

(40b) 
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(41a) 


The  X*  and  y-components  of  and  are,  in  turn,  given  by 

“nx  = ^ ~ ^ + y sin  ^ cos  ^ 

any  = 3n  sin  0 = X sin  0 cos  ^ + y sin^  0 (41b) 

= -aj  sin  0 = X sin^  4>-y  sin  0 cos  0 (4lc) 

a,y  = cos  ^ = - X sin  ^ cos  0 + y cos^  0 (41d) 

The  inertia  force  for  the  cable  segment  in  the  spatial  coordinates  x and  y takes  the  fortn 

F,c  =I'Pfio-*-“^A£n)a„^+M«najjt 

+ l(pCn +apAfin)any +M6oa,yl  T (42) 

where  i , j are  unit  vectors  in  the  x-  and  y-directions,  respectively. 

If  one-half  of  the  inertia  force  of  the  cable  segments  above  and  below  the  node  are 
summed  and  the  possible  presence  of  an  intermediate  body  at  the  node  is  taken  into  account. 
Equations  (39)  to  (42)  yield  the  following  equation  for  the  total  inert'a  force  F,  at  the  node 

= 1^ + cos*^,  + cos^^+Mgvxj  X 

+ y ; sin  cos  0^  sin  ^ cos  yj  i 

L \ ^ " 2 s«n  (>j,  cos  x 

+ + o‘p\  <>t  * ^p\  ^ob  ^ ^Bvy)  y]  T 

= (J^x  + Ky)T  + (Kx  + Jyy)7  (43) 

where  subscripts  a and  b respectively  deno:-  the  cable  segments  above  and  below  the  node 
and  Mgyj^  and  M^vy  respectively  the  virtual  i.'ass,  the  mass  plus  the  added  mass,  of  the 
intermediate  body  for  motions  in  the  x-  and  y-directions. 

Terraon  Force 

For  an  extensible  cable  segment,  the  tension  force  T depends  on  the  strain  e and  the 
strain  rate  e 
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1 = T(t.t) 


(44) 


where 


'T"-' 


(45) 


^ " dt 


'it., 


(46) 


Here  t is  the  stretched  length  of  the  cable  segment  and  the  Mihscript><  V:  and  u res|X-cti\ely 
refer  to  the  lower  and  upper  end.  of  the  cable  segment.  I 

Several  relatioaships  have  been  proposed  in  the  literature  for  th#  form  of  this  dependence. 
In  the  present  study,  the  following  relatjvely  simple  yet  general  relationship  is  used 

! 


5 = 


(47) 


where  is  the  internal  damping  A>cfltcient.  Note  that  the  ahovof dynamic  tension  differs 
from  the  static  tension.  Hquation  (X)  by  the  addition  of  the  lineilr  internal  structural  damping 
term  C,  e.  I his  model,  commonly  refera'd  to  as  the  Voigt  inodtl.  has  been  used  in  previous 
cable  studies,  for  example,  by  Huffman’*'  and  by  Gwiler  and^Cuira  Goeller  and  Laura 
show  experimentally  that  this  mode!  adequately  describes  infernal  damping  for  nylon  ropes 
except  for  frequencies  significantly  higher  than  the  resimance  frequency  of  the  cable.  Other 
more  complex  forms  for  the  internal  structural  damping,  such  as  those  proposed  by  Reid.  ’' 
can  !X’  conveniently  incorporated  into  the  program. 

Since  the  c.ible  segments  abo*  e and  K'l^jj^ktlle  node  aet  on  it.  the  lollowing  two 
expressions  are  obtained  for  and  ( rcspeelively  the  and  y -components  of  the  resultant 
fotecs  due  to  eabk  tension). 


■*v  = *ay  " ^hy  = "‘a  ^ ^ 

/ 

■’"lli.llman.  H.R...  "The  Itynaniical  Behavior  of  hxtensible  fable  in  a Unsfiinii  Flow  Field  (.4n  Invcsligaiion 
ot  ilie  lowed  Vehicle  Piohlcni)."  Pti.ti.  Thesis.  Purdue  llniversits  (Jan  P)6*)) 

ellci.  J.F  ai.d  P..A.  l.aura.  ■■Analytical  and  Fxperimenial  Simly  o'"  the  Dynamic  Response  of  fable 
Systc"'.-."  The  falholic  fnivcrsity  ol  America,  Deparlmcnl  o!  Mechanical  1 ngmeering.  ‘''hcmis  Program 
R.ro.l  7()-.t  t.Api  1^70). 

’•Rcia.  R O . ■■Dynamics  ol  l>ecp-Sea  Mo  -(ing  lanes,”  Texas  A&.M  I'niverse'. . Department  ol  Oceano- 
;.iaph>.  A/ M Protect  .’04.  Reicience  nX-IlF  (Jul  l*)68). 
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Drag  Forces 


As  shown  previously  in  Equations  (6)  and  (7),  the  static  normal  and  tangeniul  c.iblc 
drags  are  taken  to  be  respectively  proportional  to  the  squares  of  the  normal  and  tangential 
components  of  the  current  velocity  c.  In  the  dynamic  case,  the  fluid  velocity  relative  to  the 
cable  must  include  the  ocean  wave  particle  velocities  and  y^  as  well  as  the  velocities  x 
and  y of  th<*  cable.  The  resultant  expressions  and  v^^  for  the  relative  velocities  normal 
and  tangential  to  the  cable  then  take  the  form 

Vjji  = (c  + x^  - X)  cos  0 + (y^  - y)  sin  0 = x^  cos  0 + sin  0 l49al 

v^j  - -(c  + x^  - x)  sin  0 + (y^  - y)  cos  0 = - x^  sin  0 + y,  cos  0 (4%) 

where  x^  = c + - x and  y^  = - y. 

If  one-half  of  the  drag  forces  acting  on  the  cable  segments  above  and  below  the  node  are 
summed,  the  following  equations  arc  obtained  for  the  resultant  forces  and 

Djy  acting  on  the  rode 

^x  = -TDuSin0,--^D,^sin0j, 

where 

^n(a.b)  ~ D(a.b)  ‘^(a.b)  ^Ot.i.b)  'mta.bl  ^ 'rnla.b)  ^ 

^l(a.b)  ~ ~2  T(a.b)  ‘^(a.bl  ^«i<a,b)  'rtla.b) ' 'rlla.bi  ^ 

'rn(a.b)  = \ <>(a.b)  * Vbl 

't;(a.bl  = - <*(a.b)  ^a.bi 


(.‘iOal 
I.^Obl 
l.'^OcI 
( 50d» 
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INTERMEDIATE  BODIES 


The  magnitude  and  direction  of  the  added  inertia  and  drag  forces  for  an  arbitrary  body 
depend,  in  general,  in  a complex  manner  on  body  shape  and  orientation  relative  to  the  flow. 
The  modeling  of  intermediate  bodies  in  previous  cable  studies  has  ranged  from  the  very  simple 
to  the  very  complex.  At  one  extreme,  the  presence  of  intermediate  bodies  has  been  neglected 
altogether,  leading  to  a cable-only  system.  At  the  other  extreme,  some  studies  have  paid  very 
careful  attention  to  a particular  body,  often  the  lower  body  of  the  system,  and  approximated 
the  rest  of  the  cable  system  in  a simple  manner. 

After  careful  review  of  previous  studies,  formulations  were  selected  for  the  present  study 
which  although  relatively  simple,  can  model  most  bodies  of  interest  for  sonobuoy  systems. 
They  are  also  applicable  to  other  cable  systems  where  the  bodies  are  relatively  small  and/or 
conform  to  the  shape  limitations  given  below. 

CONSTANT  COEFFICIENTS 
Inirtia  Forets 


The  inertia  forces  in  the  x-and  y-directions  are  expressed  as  the  virtual  mass  times  the 
acceleration  in  these  directions,  where  the  virtual  mass  is  the  sum  of  the  mass  of  the  body 
Mg  and  the  constant  infinite  fluid  added  mass.  Thus,  the  virtual  masses  and  Mg^^ 

which  appear  in  Equation  (43)  are  given  by 

(5Ia) 

Mgvy  = Mg+KypVg 

(51b) 

where  K_,  K are  respectively  the  added  mass  coefficients  for  motions  In  the  x 

* 7 

y-directions  and  Vg  is  the  reference  volume,  usually  the  volume  of  the  body. 

- and 

Drag  Forces 

Two  formulations  are  used  to  describe  the  drag  forces.  In  one.  the  drag  components 
Dgjj  and  Dgy  are  taken  to  have  the  form 

Dbx=T^^dax\\ 

(52a) 

^By  ~~2^  ^DAy  ^r 

(52b) 
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where  C'^^y  are  respectively  the  drag  areas  in  the  x-  and  y-directions  and 

Vf  = \J + >’f“  is  the  resultant  fluid  velocity  relative  to  the  body.  This  approach  has  been 
used  by  Walton  and  Polachek^^  and  is  exact  for  the  case  of  a sphere.  In  this  case,  wrhere 
^i)Ax  “ ^DAy  ~ ^DA’  resultant  drag  Dg  is  parallel  to  the  resultant  fluid  velocity  v^ 


This  approach  may  also  be  used  to  approximate  the  drag  for  other  blunt  shapes  for  which  the 
drag  areas  for  different  flow  directions  do  not  differ  greatly,  e.g  . near-cubes  or  circular 
cylinders  with  length  to  diameter  ratios  of  approximately  1 . 

In  the  second  approach.  Dgj^  and  Dgy  are  taken  as  proportional  to  the  squares  of  the 
components  of  the  relative  fluid  velocities  in  these  respective  directions 


'bx  2 * 

(53a) 

'By  "3^^DAy  yrlVr' 

(53b) 

This  is  a good  approximation  for  long  cylinders  or  thin  disks  with  axes  parallel  to  the  x-  or 
y-directions.  In  these  cases,  where  there  is  a large  difference  between  the  drag  areas,  the  drag 
in  one  direction  is  essentially  pressure  drag  whereas  the  much  smaller  drag  in  the  other 
direction  is  evsentially  due  to  fluid  friction. 

In  the  present  program,  the  choice  of  whether  to  use  Equation  (52)  or  (53)  to  compute 
the  drag  is  determined  by  the  value  of  the  ratio  Cp^j^/Cp^y.  It  was  somewhat  arbitrarily 
decided  that  in  the  range 

0 5 < Cp^,/Cp^y  < :.0  (54) 

Equations  (52a)  and  152b)  would  be  used  to  compute  Dgj^  and  Dgy.  Outside  this  range. 
Equations  (5.'a)  and  (53b)  are  used  to  compute  Dg^^  and  Dgy. 


VARIABLE  COEFFICIENTS  FOR  CIRCULAR  DISK 

When  a ln)dy  is  executing  dynamic  oscillations  such  that  it  periodically  traverses  its  own 
VISCOUS  wake,  the  added  mas>  and  drag  coefficients  are  more  conectly  expressed  as  functions 


"Walton.  T.S.  and  H.  Polachek.  'Talcubiion  of  Transient  Motion  of  Submerged  Cables,”  Mathematical 
Tables  and  Other  Aids  to  ComputatKm.  Vol.  14.  No.  69.  pp.  27  -46  (Jan  I960). 
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os'  the  dynamic  motion.  Drag  coefncients  are  usually  obtained  for  steady  flow  and  added 
mass  coefficients  are  often  computed  for  an  inviscid  fluid.  The  reviews  by  Wiegel^  and 
Holler^^  show  that  nearly  ail  of  the  measureirents  of  the  dynamic  coefficients  for  various 
bodies  have  been  conducted  for  dynamic  oscillation  in  one  direction  only  and  in  the  absence 
of  any  steady-state  fluid  velocity.  Thus,  the  results  from  these  studies  should  be  used  with 
caution  in  the  present  case  where  the  cable  system  will  generally  undergo  dynamic  motions  in 
both  directions  in  the  presence  of  a steady-state  current  profile. 

For  the  particular  case  of  a circular  disk,  which  is  commonly  used  in  sonobucy  systems 
to  damp  out  the  motions  of  the  lower  acoustic  units,  the  user  may  either  employ  the  constant 
coefficient  approach  (described  previously)  or  have  the  program  internally  compute  the 
dynamic  added  mass  and  drag  coefficients  for  the  direction  normal  to  the  disk  plane,  as  based 
on  the  experimental  relationships  given  by  Holler.^*  The  coefficients  arc  related  to  the 
dynamic  motions  as  follows^^ 


C 


DAn 


4 y 


(55a) 


(S5b) 


= y/jM 


for  0.077  < a < 3.84 
= 1 .96  for  0 > 3.84 


= y/O^  = 0.278  for  B < 0.077 


where 


n = direction  normal  to  the  disk  plane,  either  x or  y 
= added  mass 

h.  ii  = the  velocity  and  acceleration,  respectively,  of  the  disk  in  this  direction, 
either  (x.x)  or  (y.y) 


The  above  formulas  show  how  the  drag  area  and  the  aJdcd  mass  vary  with  0.  which  is  a 
measure  of  the  ratio  of  relative  magnitudes  of  the  velocity  and  acceleration.  At  low  values  of 
0,  where  the  accelera.ion  is  much  higher  than  the  velocity  (e.g..  during  the  initial  instants  of 


^Wiegel.  R.L,  "Oceanographical  Engineering.”  Prentice  liall,  Inc..  Englewood  Giffs,  NJ..  (1964), 
Oiapter  II. 

^■'Holler.  R.A.,  “Hydrodynamic  Effects  of  Hartitonic  Acceleration.”  Naval  An  Development  Center 
Rejxirt  AE-”i:0(Jan  1 972). 
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a body  starting  from  rest),  the  formulas  show  that  has  a high  value  and  that  the  added 
mass  is  given  by  = pdg^/3,  the  potential  flow  result.  At  the  other  extreme  of  high  values 
of  p.  where  the  velocity  is  much  higher  than  the  acceleration,  the  added  mass  has  a high  value 
and  the  drag  area  is  given  by  = 1.12»dg /4,  the  steady-state  value. 

The  added  mass  and  drag  in  the  direction  tangent  to  the  disk  plane,  which  are  much 
smaller,  are  computed  by  the  constant  coefficient  approach  outlined  previously. 

DESCRIPTION  OF  COMPUTER  PROGRAM 

Program  CABUOY  consists  of  a main  program  and  six  subroutines. 

MAIN  PROGRAM 

The  main  program  accepts  input  data  for  the  cable  system,  surface  waves,  current  profile, 
and  the  initial  conditions  for  the  dynamic  calculations.  If  a surface  buoy  is  present,  input 
data  are  read  in  by  Subroutine  BUOY,  described  below.  Data  may  be  entered  in  either 
English  or  metric  units.  A detailed  description  of  input  instructions  is  given  later. 

The  program  is  currently  written  to  accept  up  to  50  cable  segments  and  4‘)  intermediate 
bodies.  This  number  can  be  conveniently  increased  by  changing  a few  DIMENSION  and 
COMMON  statements,  but  it  should  be  noted  that  dynamic  calculations  for  more  than  50 
nodes  arc  likely  to  require  prohibitively  large  amounts  of  computer  time. 

The  main  program  pnnts  out  the  input  data  and  then  calls  on  various  subroutines  to 
;alculatc  the  ocean  wave  spectrum,  certain  constants  for  the  surface  buoy,  the  steady-state 
configuration  of  the  cable  system,  and  finally  the  dynamic  motions  of  the  system  at  prescribed 
time  intervals.  The  output  from  the  steady-state  and  dynamic  calculations  are  also  printed  by 
this  program. 

SUBROUTINE  STAT 

This  subroutine  defines  the  five  steady-state  differential  equations.  ( I ) to  (5),  for  the  cable. 

SUBROUTINE  DYNA 

I his  subroutine  defines  the  dynamic  differential  equations.  (36)  and  (37).  for  each  of  the 
M nodes.  For  cases  where  a buoy  is  present,  this  subroutine  also  defines  the  three  diftcrential 
equations  of  motion  for  the  surface  buoy,  namely  (23)  to  (25). 
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SUBROUTINE  CUR 


This  Mibroutine  furnishes  the  steady-state  current  profile  relative  to  the  cable  system. 

For  a free-fioating  cable  system,  this  would  be  the  actual  current  profile  minus  the  drift 
velocity  of  the  cable  system.  For  a given  value  of  the  vertical  distance  y,  this  subroutine 
linearly  interpolates  between  the  input  velocities  which  are  read  in  as  a function  of  y.  For 
cases  where  the  given  value  of  y is  greater  (less)  than  the  largest  (smallest)  value  of  y which 
IS  read  in.  the  subroutine  takes  the  velocity  to  be  the  value  at  the  largest  (smallest)  algebraic 
value  of  y which  is  read  in. 

SUBROUTINE  SPECT 

This  subroutine  is  employed  when  the  user  wishes  the  program  to  internally  generate  the 
surface  wave  components.  In  this  case,  the  subroutine  defines  the  amplitudes  of  the  surface 
wave  components  by  using  the  Pierson-Moskowitz  energy  sea  spectrum.  Equation  (18). 
Provision  is  left  at  the  end  of  the  subroutine  for  implementing  other  forms  for  the  sea 
spectrum. 

SUBROUTINE  BuOY 

This  subroutine  is  used  when  a surface  buoy  is  present.  After  reading  input  data  for  the 
surface  buoy,  the  subroutin*’  calculates  the  various  buoy  geometrical  and  added  inertia 
coefficients  which  appear  in  Equations  \ 2~f)  to  (35).  It  concludes  by  calculating  the  steady- 
state  pitch  angle  of  the  buoy,  in  the  absence  of  any  dynamic  excitation  due  to  surface  waves. 

SUBROUTINE  ITERA 

1 his  subroutine  is  used  for  boundary  value  cases  when  iteration  schemes  are  required  to 
obtain  the  steady-state  configuration  of  the  cable  system.  It  contains  iteration  schemes  which 
are  applicable  for  free-floating  cable  systems  and  a cable  of  given  length  moored  in  a given 
ocean  depth.  Provision  is  made  at  the  end  of  the  subroutine  for  implementing  iteration 
schemes  for  other  applications. 

SUBROUTINE  KUTMER 

Ih’.s  subroutine  uses  the  Kutta-Merson  method  to  numerically  integrate  the  steady-state 
difl'erential  equations  defined  in  Subroutine  ST AT  and  the  dynamic  differential  equations  of 
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motion  defined  in  Subroutine  DYNA.  The  subroutine  automatically  reduces  the  integration 
step  size  until  specified  error  criteria  arc  met. 


INPUT  INSTRUCTIONS 


READ  STATEMENTS 


Input  data  are  entered  into  the  program  by  means  of  the  following  READ  statements 
contained  in  Program  MAIN  and  Subroutine  BUOY.  These  statements  are  given  numbers 
simply  for  identification  purposes. 

MAIN  Program 


READ  (5,1) 

NCASES 

Card  1 

DO  1000  MC 

-1.  NCASES 

READ  (5.301)  TITLE 

Card  2 

READ  (5.1) 

NSM.  NSW.  NCAB,  NCUR,  ITER.  MTRC 

Card  3 

READ  (5,2) 

(FSM(K),  K=1,NSM) 

Card  4 

READ  (5.2) 

(AXSM(K).  K=I.NSM) 

Card  5 

READ  (5.2) 

(.AYSM(K).  K=1,NSM) 

Card  b 

READ  (5,2) 

(FIDSM(K).  K=1.NSM) 

Card  7 

READ  (5.2) 

(ASW(K).  K=  l.NSW) 

Card  8 

READ  (5.2) 

(FRSW(K).  K=I.NSW) 

Card 

READ  (5.2) 

(FIDSW(K).K=  l.NSW) 

Card  10 

READ  (5,2) 

RHO.  SUBM.  TWX,  TIY,  CDASX.  AMC.  AFAC.  TMIN 

Card  1 1 

READ  (5,2) 

TlNVl , DTI . TOTT,  DT2,  DIR.  TBH,  TBYMX 

Card  1 2 

READ  (5,3) 

(FLC(K).  K=  I,  NCAB) 

Card  13 

READ  (5.2) 

(DCKK).  K=1,NCAB) 

Card  14 

READ  (5.2) 

(CDN(K).  K=1.NCAB) 

Card  1 5 

READ  (5,2) 

(CDT(K).  K=  l.NCAB) 

Card  16 

READ  (5,2) 

(WC(K).  K=  l.NCAB) 

Card  1 7 

READ  (5,4) 

(CM(K).  K=  1..NCAB) 

Card  1 8 

READ  (5.3) 

(TREF(K).K=  l.NCAB) 

Card  1 9 

READ  (5.5) 

(Cl  (K).K=  l.NCAB) 

Card  20 

READ  (5.2) 

(C2(K).  K = 1,  NCAB) 

Card  21 

READ  (5.2) 

(C1NT(K),K=  l.NCAB) 

Card  22 

READ  (5.2) 

(WBD(K),K=  l.NCAB) 

Card  23 

READ  (5. 2) 

(CDABX(K).  K=  l.NCAB) 

Card  24 

READ  (5.2) 

(CDABY(K).K=  l.NCAB) 

Card  25 

READ  (5,2) 

(XMBV(K).  K=  l.NC.AB) 

Card  2h 

READ  (5.2) 

(YMBV(K).  K=  1.  NCAB) 

Card  27 

READ  (5.3) 

(YY(I).  1=1.  NCUR) 

Card  28 

READ  (5.3) 

(CCK(I).  I=1.NCUR) 

Card  29 

READ  (5.2) 

(PHID(I).  1 = l.NCAB) 

Card  30 

READ(5.3)  (TENKI).  I=I,NCAB) 

READ (5.2)  (XPKI),  1=1,NCAB) 

READ  (5.2)  (YPI(I).I  = 1.NCAB) 

1000  CONTINUE 

The  corresponding  FORMAT  statements  are: 

1 FORMAT  (2413> 

2 FORMAT  (8F10.4) 

3 FORMAT  (8F102) 

4 FORMAT  (8FI0.6> 

5 FORMAT  (8F10.01 
301  FORMAT  (20A4) 

Subroutine  BUOY 

READ  (5.1)  CDASY,  WAS,  RWY.  RTX.  RTY.  YCG.  BIN  (,ard34 

READ(5.n  XSI.  ZETI.  SYDl.  XPSl.  ZTPI.  SYPDl  (ard  35 

The  corresponding  FORMAT  statement  is: 

I FORMAT  (8F10.4) 


Card  3 1 
Caid  32 
Card  33 


DEFINITION  OF  INPUT  VARIABLES  FOR  MAIN  PROGRAM 


NCASES 

TITLE 

NSM* 

NSW^ 

NCAB 

NCUR 

MTRC 

ITER 

FSM(K)" 

AXSM(K)- 

AYSM(K)- 

FIDSM(K)- 

ASW(K)^ 


Number  of  cases.  NCASES  > I 
Title 

Number  of  surface  motion  components.  1 < NSM  < 20 

Number  of  surface  vave  components.  1 < NSW  < 20 

Number  of  cable  segments,  2 < NCAB  < 50 

Number  of  current  profile  points.  2 < NCUR  < 10 

MTRC  < 0 if  input  data  are  entered  in  English  units;  MTRC  > 1 it  input 
data  arc  entered  in  metric  units 

Iteration  index 

NSM 

XsM  = ^ AXSM(Kfcos(-2ir*  FSM(Kt*t  + FIDSMtKftr  180.1 
k =1 

NSM 

VsM  ^ 1C  -AYSMtKfsin  (- 2ir*  FSM(Kl*t  + FIDSMtKfir  !80.i 

k = I 
NSW 

^sw  =C  A,SW<K)*cos(-2jr*  FR.SW(Kft  + FlDSW(Kl*rr  180.) 
k = 1 


33 


FRSW(K)^ 

F!DSW(K)^ 

RHO 

SUBM< 

TWX^ 

TIY 

CDASX^ 

AMC 

AFAC 

TMIN 

TIN  VI 

DTI 

TOTT 

DT2 

DIR 

TBH 

TBYMX 

FLC(K) 

DCKK) 

CDN(K) 

CDT(K) 

WOK) 

CM(K) 

TREF(K) 

CI(K)\.C2(K). 

CINT(K) 

WBD(K) 

CDABX(K)*. 

CDABY<K)‘’ 

XMBV(K)^ 

YMBV(K)^ 


NSH 

ysw=5D  -ASW(K)*sin  (- 2ir‘ FRSW(K)*t  + FIDSW(K)*7r/l 80.) 
k = \ 

Fluid  density  in  slugs/feet^  ( kilograms/meters^ ) 

Submergence  of  top  point  of  cable  below  free  surface  in  feet  (meters) 

Horizontal  force  acting  at  top  of  cable  in  pounds  (newtons) 

Vertical  component  of  tension  at  top  of  cable  in  pounds  (newtons) 

Drag  area  of  surface  buoy  perpendicular  to  the  x-axis  in  feet^  (meters^) 

Added  mass  coefficient  of  cable;  AMC  = 1.0  :or  round  cable 

Cross-sectional  area  of  cable  = AFAC*jrd^  /4;  AFAC  = 1.0  for  round  cable 

Minimum  algebraic  tension  which  can  oe  supported  b>  c ible  in  pounds 
(newtons) 

Initial  time  interval  in  seconds  for  dynamic  calculations 

Time  step  in  seconds  for  which  printout  is  desired  for  0 < t < TINVI 

Total  time  in  seconds  for  which  dynamic  calculations  are  desired 

Time  step  in  seconds  for  which  printout  is  desired  for  TINV 1 < t < TOTT 

DIR  < 0.  if  initial  conditions  are  prescribed  at  the  bottom  (towing  cable 
case);  otherwise  DIR  > 0 

Applied  force  in  pounds  (newtons)  on  lower  weight,  body  NCAB-i.  in  x- 
direction 

Maximum  absolute  value  in  pounds  (newtons)  of  tension  in  cable  just  below 
buoy;  for  buoy-cable  system,  set  TBYMX  equal  to  a large  number,  say.  90999 

Length  of  Kth  cable  segment  in  feet  (meters) 

Diameter  of  Kth  cable  segment  in  inches  (centimeters) 

Normal  drag  coefficient  of  Kth  cable  segment 

Tangential  drag  coefficient  of  Kth  cable  segment 

Weight  in  fluid  in  pounds/foot  (newtons/meter)  of  Kth  cable  segment  at  the 
reference  cable  tension 

Mass  of  Kth  cable  segment  in  slugs/foot  (ki'ograms ’meter)  at  the  referer  'c 
cable  tension 

Reference  tension  in  pounds  (newtons)  of  Kth  cable  segment 

Tension  = TREF(K)  + Cl  (K)*e^^‘*^*  + ClNT(K)*e;  for  linearly  elastic 
material.  Cl  (K)  = AE  and  C2(K)  = 1 

Weight  in  fluid  of  Kth  body  in  pounds  (newtons) 

Drag  area  of  Kth  body  in  feet^  (meters^ ) for  flow  in  (x.y)  directions 

Virtual  mass  (mass  -Padded  mass)  in  slugs  (kilograms)  of  Kth  body  in  (x.y) 
directions 
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YY(I)’ 

CCK(I) 

PH1D(I)« 

TENKD* 

XPKI) 

YPKl) 


Value  of  y in  feet  (meters) 

Value  of  current  in  knots  (meters/second)  at  y = YY(I) 

Initial  value  of  4>  of  1th  cable  segment  in  degrees 
Initial  value  of  tension  of  Ith  cable  segment  in  pounds  (newtons) 
Initial  value  of  x of  Ith  node  in  feet/second  ( metcrs/second ) 
Initial  value  of  y of  Ith  node  in  feet/second  (meters 'second) 


DEFINITION  OF  INPUT  VARIABLES  FOR  SURFACE  BUOY 


CDASY 

WAS 

RWY 

RTX.  RTY 
YCXJ 

BIN 

XSI.  ZETl. 
SYDI’ 


Drag  area  for  y-direction  in  feet*  (meters* ) 

Weight  in  air  in  pounds  (newtons) 

Vertical  distance  of  wind  loading  center  of  pressua*  from  buoy  center  of 
gravity  Y(X  in  feet  (meters) 

(x.y)  distance  of  cable  attachment  point  from  YCG  in  feet  (meters) 

Submergence  of  center  of  gravity  below  the  free  surface  under  the  action  of 
its  own  weight  in  air  WAS  and  the  vertical  component  of  the  steady-state 
tension  (-T1Y)  in  feet  (meters) 

Moment  of  inertia  in  air  about  YCXj  in  slug  feet*  (kilogram  meters* ) 

initial  values  of  (x.  ti>)  in  (feet.  feet,  degrees)) meters,  meters,  degrees), 
where  f is  the  vertical  displacement  of  the  a'nter  of  gravity  from  its 
equilibrium  value  YCXI 


XPSI,  ZTPI. 
SYPDl 


Initial  values  of  (x.  j.  i,>)  in  (fect/sccond.  feet 'second,  degrees 'second ) 
(meters/second,  meters/sccond.  degrees 'second) 


EXPLANATORY  NOTES 

1.  ITER  = 0.  no  iteration  (prescribed  initial  steady-state  conditions) 

1 . free-floating  cable  system 

2.  moored  cable  with  given  length  in  given  depth 

3.  iteration  schenw  to  be  programmed  by  user 

2.  For  1000.  < FSM(l)  < 2000..  the  program  makes  the  presenbed  surface  motion 
components  equal  to  the  surface  wave  components  by  setting  AXSMtK)  = AYSM(K) 

= ASW(K).  FSM(K)  = FRSW(K).  and  FIDSM(K)  = FIDSW(K)  for  K = I to  K = NSM; 
the  program  automatically  sets  NSM  = NSW. 

Fo*  2(XX).  < FSM(  I ) < 30(X)..  the  program  accepts  input  data  for  a spar  buoy  and 
considers  AXSM(K)  to  be  the  cross-sectional  area  of  the  buoy  in  feet*  (meters* ) at  depth 
AYSM(K)  feet  (meters*  below  the  tree  surface.  .AYSM(  1 ) = 0.  and  .-NYSMlNSM)  = total 
draft  under  the  combined  action  of  buoy  weight  in  air  and  the  vertical  component  of  the 
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steady-state  tension.  NSM  should  be  an  odd  number.  The  input  values  for  F'DSMtK)  may 
take  on  any  values  such  as,  say.  0. 

For  FSM(  1 ) > 3000..  the  program  accepts  input  data  for  a spheroidal  buoy  and  considers 
AXSM(  I Mo  be  the  radius  of  the  buoy  cross  section  at  the  free  surface  and  AYSMt  1 ) to  be 
the  total  draft.  The  rest  of  the  input  values  of  AXSM(K)  and  AYSM(K)  as  well  as  all  of  the 
FIDSM(K)  may  tak**  on  any  values,  e.g..  0. 

3.  For  .ASW(  I ) > 1000..  the  program  computes  the  amplitude  of  the  .ASW  surface  wave 
co  mponents  by  using  the  Pierson-Moskowitz  sea  spectrum.  In  these  cases,  the  program 
considers  the  significant  wave  height  in  feet  (meters!  to  be  (ASW(  1 ) - lOOO.l  and  FRSW( ! » 
and  FRSW(21  to  respectively  be  the  lower  and  upper  frequencies  of  the  spectrum  in  cycles 
per  second.  The  program  internally  genemtes  the  phases  of  the  wave  components  by  consider- 
ing them  to  be  uniformly  separated  by  360,/NSW  degrees.  The  phase  of  the  lowest  frequency 
component,  in  degrees,  is  taken  to  be  the  input  value  of  FIDSW(  1 ). 

4.  For  the  case  of  a surfat'c  buoy  (FSM(  1 » > 2000.).  the  program  calculates  the  drag  acting 
on  the  surface  buoy  due  to  the  ocean  current  by  taking  the  value  of  the  ocean  current  SUBM 
feet  (meters)  below  the  free  surface.  Thus.  0 ^ SUBM  < total  draft. 

The  total  horizontal  forex*  at  the  top  point  of  the  cable  TIX  = TWX 
+ ( 1 lip*  CDASX*  C’C'F  (SUBM)*  .'BS  (CCF  (SUBM)).  In  cases  where  there  is  no  surface 
buoy  (i.e..  prescribed  surfaa*  motion),  TWX  and/or  CDASX  may  be  set  equal  to  zero.  For 
cases  of  a surface  buoy.  TWX  represents  the  wind  loading  on  the  buoy  in  pounds  (newtons). 

5.  For  free-floating  and  towing  cables  where  the  last  (K  = NCAB)  cable  connecting  the  lower 
weight  to  the  ocean  lx)ltom  is  fictitious,  read  in  a value  for  Cl  (NCAB)  less  than  0.0001 
pounds  (0.0004  newtons).  In  these  cases,  the  program  sets  IX'I(NCAB)  = CDN(NCAB) 

= CDT(NCAB)  = WC  (NCAB)  = CM(NCAB)  = CINT(N(  AB)  = 0.  FLC(NCAB) 

= 2 ' FLC(NCAB-  I ).  and  ( 2(NCAB)  = I. 

(>.  If  CD.-\BX(K)  IS  negative,  the  program  considers  the  body  to  be  a circular  disk  with 
plane  perpendicular  to  the  x-a\is  and  calculates  drag  and  ailded  mass  forces  by  using  the 
formulation  given  in  liquation  (5.*').  In  these  cases.  Cl)ABX(K)  is  the  negative  of  the  actual 
drag  area  and  XMBV'(K)  is  the  mass  (not  the  virtual  mass)  of  the  disk.  In  these  cases. 
CI)ABY(K)  and  YMBV(K)  should  l>e  ptvdtivc  and  retain  the  definitions  given  previousK. 
Similar  remarks  apply  if  Cl)ABY(K)  is  read  in  as  a negative  niimlwr  except  that  the  plane  of 
tlie  disk  is  now  perpendicular  to  the  y-axis. 

7.  When  IThR  = 2.  the  porgram  takes  YY(NCUR)  to  K*  tlie  ocean  depth. 

X.  F'or  1 PMID(  1 ) I > 3(>0..  the  program  takes  the  initial  values  «>l  the  angle  .*"•!  tension  ol 
e.ich  cable  segment  to  correspond  to  their  a'spective  steady-state  values  at  the  midpoint  ol 


3(> 


...... 


each  segment.  These  steady-state  values  have  been  previously  calculated  by  the  pregram. 

This  approach  will  minimize  transient  dynamic  effects.  In  these  cases,  input  values  for  the 
remaining  PHID(K)  as  well  as  all  of  the  TENI(K)  may  be  arbitrary,  e.g  0. 

9.  For  SYDI  > 360.,  the  program  sets  the  initial  value  for  buoy  inclination  ^ equal  to  the 
steady-state  value  of  4>,  which  has  previously  been  calculated  by  the  program.  This  will  tend 
to  minimize  transient  dynamic  motions  of  the  surface  body. 

PROGRAM  STORAGE  AND  TIME  REQUIREMENTS 

On  the  CDC  6700  currently  in  use  at  the  Onter,  the  program  requires  t memory  of 
approximately  47,200  octal  words  to  load  and  33,700  octal  words  to  execute.  Compilation 
time  is  approximately  23  seconds.  Execution  time  for  a particular  case  depends  on  a large 
number  of  factors,  the  most  important  of  which  include  the  tehsion-s.train  relation  of  the 
cable,  the  number  of  cable  nodes,  the  frequencies  of  the  exciting  surface  viraVes  and  the 
prescribed  surface  motion  (if  any),  and  the  amount  of  time  over  which  the  dynamic  motions 
are  desired.  Table  2 shows  the  computer  execution  time  ET  and  cost  for  all  of  the  sample 
problems  presented  in  the  following  chapter.  A computer  priority  (CP)  of  P2  indicates  over- 
night priority  whereas  P3  is  the  standard  daytime  priority  at  the  DTNSRDC  Computer  Center. 
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TABLE  2 - C0M1*UTER  EXECUTION  TIMES  AND  COST 
FOR  ALL  THE  SAMPLE  PROBLEMS 


Cl  Total 


lb 

ET 

Cost 

Prob. 

NCAB-1 

(4.45N1 

C2 

sec 

$ 

1A 

4 

2.4  X 10’ 

1.0 

61.7 

P2 

9.93 

IB 

4 

2.4  X 10^ 

1.0 

57.8 

P2 

9.56 

1C 

4 

2.4  X 10^ 

1.0 

62.4 

P2 

9.99 

ID 

4 

2.4  X 10^ 

1.0 

85.0 

P2 

12.00 

IE 

4 

2.4  X 10® 

1.0 

250.0* 

P2 

28.00 

IF 

4 

2.4  X 10^ 

0.5 

1000.0* 

P2 

100.00 

1G 

4 

2.4  X 10^ 

2.0 

55.2 

P2 

8.93 

2A 

1 

2.4  X 10^ 

1.0 

9.3 

P2 

5.15 

2B 

2 

2.4  X 10^ 

1.0 

25.1 

P2 

6.60 

2C 

4 

2.4  X 10^ 

1.0 

62.4 

P2 

9.99 

20 

8 

2.4  X 10^ 

1.0 

180.0* 

P2 

20.00 

2E 

15 

2.4  X 10^ 

1.0 

600.0* 

P2 

60.00 

3A 

4 

2.0  X 10’ 

1.0 

90.0** *** **** 

P3 

15.00 

3B 

4 

2.0  X 10’ 

1.0 

100.0** 

P3 

16.00 

3C 

4 

2.0  X 10’ 

1.0 

100.0** 

P3 

16.00 

4A 

4 

82.6 

P3 

14.27 

4B 

4 

• ••• 

#♦«« 

74.2 

P2 

11.10 

4C*‘* 

4 

• ••• 

• ••• 

100.0* 

P3 

16.00 

* Extrapolated  to  50  sec  of  dynamic  motion 

**  Extrapolated  to  20  sec  of  deployment  time 

***  Six  surface  wave  components 

****  See  Figure  11 


SAMPLE  PROBLEMS 

Input  cards  are  listed  for  the  four  sample  problems  presented  to  illustrate  use  of  the 
program.  Representative  portions  of  the  prtsgram  output  arc  listed  for  one  of  the  cases  of 
Problem  1.  but  x^nly  sttme  final  results  arc  shown  for  the  other  three  problems. 


PROBLEM  1 - UNIFORM  CABLE  WITH  VARIOUS  TENSION-STRAIN  RELATIONS 

Problem:  Compute  the  dynamic  behavior  of  a cable  suspended  from  an  ocean  platform  in 
the  presence  of  a uniform  current  of  1 knot  in  the  +x-direction  for  different  values  of  the 
elastic  constants  Cj  and  C,  appearing  in  Equation  (47» 
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Case 

Cl 

lb 

(4.4.SN) 

C2 

A 

2.4  X 10 

1.0 

B 

2.4  X 10- 

1.0 

C 

2.4  X lO’ 

1.0 

D 

2.4  X 10-* 

1.0 

E 

2.4  X 10* 

1.0 

F 

2.4.x  10* 

O.-S 

G 

2.4  X 10* 

2.0 

The  fixed  surface  motion,  surface  wave,  cable,  and  lower  body  parameters  are  as  lollows. 


Surface  motion  = surface  wave: 
number  of  components 
frequency 

amplitude  in  x-  and  y-directions 
phase  angle 

Cable: 

length 
diameter 

normal  drag  coefficient 
tangential  drag  coefficient 
weight  in  fluid 
mass 

reference  tension 
internal  damping  coefficient 
Lower  weight: 

weight  in  fluid 
drag  area  in  x-direction 
drag  area  in  y-direction 
virtual  mass  in  x-direction 
virtual  mass  in  y-direction 
Fluid  density  : 

Represent  the  cable  by 
below  the  lower  weight, 
motions,  i.e.  .SO  sec.  For  the  initial  interval 


1 

0.1  cps 

10  ft  I.^.OS  m) 

0 deg 

1000  ft  (.^05  m) 

0.2  m.  (0.50X  cm» 

M 

0.0 :. 

0.01  lb  ft  (0  14(>  N ml 
0.001  slugs  ft  10.0478  kg  m) 

25  lb  (1 11. 2 Nt 

0 

20  lb  (8‘>  N) 

0.3  ft-  t0.027‘>  m-) 

0. 3  ft-  I0.027O  m-t 
1 .0  slugs  1 14.()  Kgt 

1 .0  slugs  ( I4.f>  kgl 

1. ‘)4  slugs  ft’  ( lOOO.O  kg  m’ » 

not  including  the  fictitious  cable  segment 

c\  cles  «if  the  surface 
if  10  sec.  print  out  the  transient  motions  every 


four  equal  segments. 

Dvnamic  motions  are  desired  lor  a total  ot  live 


3'J 


0.25  sec.  For  tl'.e  final  40  sec,  increase  the  printout  interval  to  1.0  sec.  Let  the  cable  start 
from  rest  with  the  initial  angle  and  tension  of  each  segment  equal  to  the  steady-state  values. 
Solution:  The  data  cards  for  this  problem  are  listed  in  Table  3.  The  cards  which  are  the  same 
for  all  the  cases  are  listed  at  the  top  of  the  tabl^.  The  cards  for  the  title  and  the  elastic 
constants  Cl  and  C2.  which  differ  for  each  case,  are  listed  at  the  bottom  of  the  table.  The 
symbol  b denotes  a blank  in  this  and  subsequent  tables  which  list  data  cards  for  the  sample 
problems.  Also,  Column  1,  11,  21.  31,  41,  51,  61,  and  71  have  been  indicated  since  most  of 
the  data  start  in  these  columns. 

TABLE  3 - INPUT  DATA  FOR  SAMPLE  PROBLEM  1 


Card  1 
Card  3 
Card  4 
Card  5 
Card  6 
Card  7 
Card  8 
Card  9 
Card  10 
Card  11 
Card  12 
Card  13 
Card  14 
Card  15 
Card  16 
Card  17 
Card  18 
Card  19 
Card  22 
Care'  23 
Card  24 
Card  25 
Card  26 
Card  27 
Card  28 
Card  29 


1 

11 

bbl 

bb1bb1bb5bb2bb0 

0.1 

10.0 

10.0 

0. 

10.0 

0.1 

0. 

1.94 

0. 

10. 

0.25 

250. 

250. 

0.2 

1.4 

1.4 

0.02 

0.02 

0.01 

0.01 

0.001 

0.001 

25. 

25. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

10000. 

1. 

1. 

TABLE  3 — (continued I 


1 

11 

21 

31 

41 

51 

61 

71 

Card  30 

999. 

0. 

0. 

0 

Card  31 

0. 

0. 

0. 

0. 

Card  32 

0. 

0 

0. 

0. 

Card  33 

0. 

0. 

0. 

0. 

Card  2 
Card  20 
Card  21 


bbb  PROBLEM  1A. 
24.  24. 

1 1. 


Cl  = 24,  C2  = 1 
24.  24. 

1.  1. 


Card  2 
Card  20 
Card  21 


bbb  PROBLEM  IB, 
240.  240. 

1.  1. 


Cl  = 240,  C2  = 1 
240.  240. 

1.  1. 


Card  2 
Card  20 
Card  21 


bbb  PROBLEM  1C. 
2400.  2400. 

1.  1. 


Cl  = 2400.  C2  = 1 
2400.  2400. 

1.  1. 


Card  2 
Card  20 
Card  21 


bbb  PROBLEM  10. 
24000.  24000. 

1.  1. 


Cl  = 24000.  C2  = 1 
24000.  24000. 

1.  1. 


Card  2 
Card  20 
Card  21 


bbb  PROBLEM  IE, 
240000.  240000. 

1.  1. 


Cl  = 240000.  C2  - 1 
240000.  240000. 

1.  1. 


I Card  2 
I Card  20 

Card  21 


bbb  PROBLEM  IF. 
2400.  2400. 

0.5  0.5 


Cl  = 2400.  C2  = 0.5 
2400.  2400. 

0.5  0.5 


Card  2 
Card  20 
Card  21 


bbb  PROBLEM  1G. 
2400.  2400. 

2.0  2.0 


Cl  = 2400.  C2  = 2 
2400.  2400. 

2.0  2.0 


Tabic  4 (he  t'lrsl  six  and  the  last  two  pages  of  tlic  a'lnpiitcr  output  lor  Problem 

II).  The  pages  which  have  been  left  out  simpK  contain  output  lor  the  dxnamic  motions  for 
intermediate  time  inicrxal'..  T.ibic  4 showN  that  the  first  page  of  the  output  lists  a i.-Me  of 
coiuersion  Irom  Lughsh  to  metric  units.  The  si'cond  page  lists  the  input  data.  l or  the 
presiMit  case  ol  a suspended  cable,  where  the  initial  conditions  are  known  at  the  lower  Nidy. 
two  integrations  are  perlormed  for  the  nIc.  ..-stale  configuration  (see  section  on  Initial  Value 
Casv.i.  I he  results  of  thes<*  two  integrations  are  given  m the  next  two  pages.  The  remainder 
ol  the  output  i^  a listing  ol  the  dynamic  displacements,  \elocities.  and  accelerations  of  the 
surt.ice  waves.  uptH’r  cable  point,  and  each  runle  at  the  presenbed  time  intervals  The  output 
also  lists  the  angle,  angular  velocity,  tension,  strain,  and  strain  rate  of  each  cable  segment. 
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INITIAL  TINE  INTENVAL*  lO.OSOO  SEC  TINE  STEM  .2509  SEC 

TOTAL  TINE*  59.0903  SEC  TINE  STEP*  1.9100  SEC 


TABLE  4 — (continued) 


• •••• 

•*9 

M n ^ 

• •••• 

m Hi*-  • M 

• •••• 
• •••• 

sst3:: 

MMMMM 

MMMMM 

3<332 

MMMMM 

MMMMM 

K M •«  M M 

M fte  M «#  M 
M M M M «• 

MMMMM 
MM  M mM 
MM  M mm 

ft.  M M M M 
ft.  4 •#  M M 
M •«  M M M 
M M M M ft* 

MMMMM 

MMMMM 

» M J»  # M 
••  9 

stSs:; 

• • • • • 
till 

^ M • 

«• 

1 

i:ss: 

4 4mk  4 

*sssS 

M M ••  M «• 
MMMMM 

mmmmm 

MMMMM 

MMMMM 

M K M 

M M M M 
Mite  MMM 

MMMMM 

MMMMM 

MMMMM 

M M?7^ 

MMMMM 
M M M M M 

MMMMM 
M M M M •« 

*1  O M M 1^ 
Ite  M M M M 
K ^ M •«  M 

III  I I 

M M 4 ••ftte 

«l  M M M 4 

1 ft  • 

«• 

1 

• 

MMMMM 

MMMMM 

M • . 

1 

4 M M M 4 

M III 

1 

M M M M M 

«#  III 

1 

M M M M tft 
M O M •«  M 
• < • • • 
««  III 

1 

ft.  • IW  M M 

•ft  1 f 1 1 
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• • • « • 
1 • 1 1 
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• • • • • 
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111! 

• • • • • 
MMMMM 

4 M « MM 

1 1 T 1 

M M M 4 M 

MMMMM 
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M M ••  M M 
• I • • • 

M M M M M 
MMMMM 
1 1 1 1 

M M 4 IW  M 

M M « ft#4 
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K • « 

•«  *t 

• • « • • 
•««•«««• 
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M M • M 

M M M M 

MMMM 
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4 M M M 
M MM  M 

• • I • 1 
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M M M ft« 

M M •«  M M 
• ft  ft  ft  • 
M M M ftte 
M <•  M ^ 

«•  ««  « ^ •* 

• « • » ^ 

^ M 

t<i  K 

M M • tfk  # 

M 

M««»M 

M MM  «•  M 

M K.  M M 9 

M M M l>»  K 

M M N M M 

M M M ••  M 
M MMMM 

M K M M M 
M MMMM 

M Ite  fw  M M 
M MMMM 
M K.  M M M 
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M MMMM 
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Figure  5 ^how^  the  tension  in  the  upper  cable  segment  for  0 < t < 20  for  Cases  1 A to 
IL,  where  Cl  increases  from  2.4  \ 10  lb  (1.07  x 10"  N)  to  2.4  \ lOMb  ( 1.07  x 10^  N).  As 
would  be  expected,  the  figure  shows  that  lension  fiuctuations  increase  with  increasing  values 
of  Cl.  At  the  highest  value  of  Cl  = 2.4  x 10^  lb  (1.07  x 10^  N).  Figure  5 shows  that  the 
tension  is  zero  at  t = 0.25  sec  (i.e..  the  cable  in  slack)  and  then  jumps  to  a value  of  89.47  |b 
(398.1  N)  at  t = 0.5  sec.  This  sudden  jump  in  tension  after  a slack  condition  is  often  referred 
to  as  “snap  loading."  At  the  other  extreme,  the  tension  shows  a pea k-to-t rough  lluctuation 
of  less  than  2.5  lb  (1 1.1  N)  at  the  lowest  value  of  Cl  = 2 4 .x  10  lb  ( 1.07  \ lo*  N). 


0 2 4 6 8 10  12  14  16  18  20 

TIME  (SEC) 

Figure  5 — Tension  in  Up|)er  Cable  Segment 
f«>r  Various  Values  of  Cl 

Fable  2 indicates  that  the  computer  execution  time  remains  at  approximateK  (>0  sec  tor 
C.ises  1 A.  IB.  and  1C  f.  r which  Cl  < 2.4  x 10-^  lb  ( 1.07  x 10'*  N).  At  higher  sallies  of  Cl. 
Cases  ID  and  IE.  execution  tune  increases  with  increasing  Cl.  F.xecution  time  also  increases 
with  decreasing  values  of  C2.  Case  IF.  The  reasons  tor  these  increases  m.iy  be  most 
conveniently  explained  by  considering  Equation  i9i  which  shows  that  tor  a given  value  ot 
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(T  - Tj, ).  t decreases  with  increasing  values  of  Cl  and  decreasing  values  ot  C2.  Small  values 
of  e.  in  turn,  mean  that  the  x-  and  y-displaccments  must  be  calculated  with  greater  precision, 
leading  to  smaller  integration  time  steps. 

PROBLEM  2 - UNIFORM  CABLE  REPRESENTED  BY  VARIOUS  NUMBERS  OF  SEGMENTS 

Problem;  Solve  Problem  1 1)  with  the  single  exception  of  representing  the  cable  by  the  following 
number  of  nodes.  NCAB-1 ; 1 . 2.  4.  8.  and  1 5.  Consider  the  nodes  to  be  equally  spaced. 
Solution;  The  data  cards  for  this  problem  arc  listed  in  Table  5.  Cards  1 and  4 to  12.  which  arc 
identical  to  those  for  Problem  1 D.  are  omitted.  Also  omitted  arc  Cards  14  to  .^3  which  are  simi- 
lar to  those  for  Problem  1 D w'ith  the  exception  that  the  number  of  entries  depends  on  NCAB. 


TABLE  5 - INPUT  D ATA  FOR  SAMPLE  PROBLEM  2 


1 11  21  31  41  51  61  71 

Card  2 

bbb  PROBLEM  2A,  NCAB  1 = 1 

Card  3 

bb1bb1bb2bb2 

Card  13 

1000. 

Card  2 

bbb  PROBLEM  2B.  NCAB-1  = 2 

Card  3 

bb1bb1bb3bb2 

Card  13 

500.  500. 

Card  2 

bbb  PROBLEM  2C.  NCAB  1 = 4 

Card  3 

bb1bb1bb5bb2 

Card  13 

250.  250.  250.  250. 

Card  2 

bbb  PROBLEM  2D.  NCAB-1  - 8 

Card  3 

bb1bb1bb9bb2 

Card  13 

125.  125.  125.  125.  125.  125.  125.  125. 

Card  13 

bbb 

Card  2 

bbb  PROBLEM  2E.  NCAB-1  = 15 

Card  3 

bb1bb1b16bb2 

Card  13 

66  667  66.667  66  667  66.667  66.667  66.667  66.667  66.667 

Card  13 

66.667  66.667  66.667  66.667  66.667  66.667  66.667 

Figure  (>  shows  how  the  steady-state  location  of  the  lower  weight  varies  with  NC.AB-I 
The  results  arc  pmir  lor  NC.AB-l  < 2.  good  for  NCAB-l  = 4,  and  converged  for  NCAB-I  > 8 
Figure  7 shmvs  the  velocities  x and  y of  the  lower  unit  at  t = 4 and  10  for  various  values 
ol  NCAB-I.  I he  results  suggest  that  a minimum  ot  four  nodes  is  required  to  obtain  maximum 
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velocities  with  accuracies  to  within  10-20  percent.  The  magnitudes  of  the  computed 
veloi,ities  generally  increase  with  increasing  number  of  nodes. 

Table  2 shows  the  manner  in  which  the  computer  execution  time  ET  increases  with  tlie 
number  of  nodes.  The  execution  times  are  approximately  proportional  to  NCAB  square  ! 

ETtt(NCAB)^  (=i(.> 


PROBLEM  3 - DEPLOYMENT  OF  SUSPENDED  CABLE  IN  CIRCULATING 
WATER  CHANNEL 


Problem:  Compute  the  deployment  of  an  initially  vertical  cable,  suspended  in  the  DTNSRIX' 
Circulating  Water  Channel,  to  its  final  steady-state  configuration  in  the  presence  of  channel 
How  speeds  of  1,  2,  and  4 knots  (0.515.  1.0.1.  2.06  m/s).  The  fixed  cable  and  lower  body 
parameters  are  as  follows: 

Cable: 


length 

diameter 

normal  drag  coefficient 
tangential  drag  coefficient 
weight  in  fluid 
mass 

reference  tension 
<^'i 

C, 

internal  damping  coefficient 
L/5wer  weight: 

weight  in  fluid 
drag  area  in  x-direction 
drag  area  in  y-direction 
virtual  mass  in  \-direction 
virtual  mass  in  y-direction 
Fluid  density: 


12  ft  (3.66  m) 

0.1 2 in.  (0.305  cm) 

1.4 

0.02 

0 

0. 0002  slugs/ft  (0.00957  kg/m) 
2 lb  (8.9  N) 

20  lb  (89.  N) 

1. 

0. 

2 lb  (8.9  N) 

0.05  ft‘  (0.00465  m-) 

0.05  ft"  (0.00465  m*) 

0.1  slugs  ( 1.46  kg) 

0.1  slugs  (1.46  kg) 

1.94  slugs/lt-  (1000.6  kg/m^  ) 


Represent  the  cable  by  four  equal  segments,  not  including  the  fictitious  cable  segment 
below  the  lower  weight.  Take  the  total  time  interval  to  be  20  sec.  For  the  initial  intcu.il  ol 
I sec.  printout  of  the  transient  cable  configuration  is  desired  every  0.1  sec.  For  the  fin.il  19 
sec.  increase  the  printout  interval  to  0.2  sec.  Let  the  system  start  from  rest  with  the  iniii.il 
tension  in  each  cable  segment  equal  to  2 lb  (8.9  N). 
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Solution:  The  data  cards  for  this  problem  are  listed  in  Table  6.  The  cards  which  are  the 
same  for  ail  the  cases  arc-  listed  at  the  top  of  the  table.  The  cards  for  the  title  and  the  current 
magnitude  which  differ  for  each  case  are  listed  at  the  bottom  of  the  table. 


TABLE  6 - INPUT  DATA  FOR  SAMPLE  PROBLEM  3 


— 

1 

11 

21 

31 

41 

^1 

61 

71 

Card  1 

bbl 

Card  3 

bb1bb1bb5bb2 

Card  4 

0.1 

Card  5 

0. 

Card  6 

0. 

Card  7 

0. 

Card  8 

0. 

Card  9 

0.1 

Card  10 

0. 

Card  11 

1.94 

0. 

0. 

-15. 

0. 

1.0 

1.0 

0. 

Card  12 

1. 

0.1 

20. 

0.2 

-1. 

0. 

99999. 

Card  13 

3. 

3. 

3. 

3. 

Card  14 

0.12 

0.12 

0.12 

0.12 

Card  15 

1.4 

1.4 

1.4 

1.4 

Card  16 

0.02 

0.02 

0.02 

0.02 

Card  17 

0. 

0. 

0. 

0. 

Card  18 

0.0002 

0.0002 

0.0002 

0.0002 

Card  19 

2. 

2. 

2. 

2. 

Card  20 

20. 

20. 

20. 

20. 

Card  21 

1.0 

1.0 

1.0 

1.0 

Card  22 

0. 

0. 

0. 

0. 

Card  23 

0. 

0. 

0. 

2. 

Card  24 

0. 

0. 

0. 

0.05 

Card  25 

0. 

0. 

0. 

0.05 

Card  26 

0. 

0. 

0. 

0.1 

Card  27 

0. 

0 

0. 

0.1 

Card  28 

0. 

100. 

Card  30 

0. 

0. 

0. 

0. 

Card  31 

2. 

2. 

2. 

2. 

Card  32 

0. 

0. 

0. 

0. 

Card  33 

0. 

0. 

0. 

0. 
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TABLE  6 — (continued) 


r 


V 

1 11  21  31  41  51  61  71 

L- 

Card  2 

bbb  PROBLEM  3A,  V = 1 )fNOT 

c 

Card  29 

1.0  1.0 

Card  2 

bbb  PROBLEM  3B,  V = 2 KNOTS 

Card  29 

2.0  2.0 

i 

Card  2 

bbb  PROBLEM  3C.  V = 4 KNOTS 

Card  29 

o 

o 

I 


k 

I Figures  X,  ‘).  and  10  show  the  configuration  of  the  cable  nodes  at  various  times  during 

^ deployment  for  currents  of  1.  2.  and  4 knots  (0.515.  l.O.L  2.0b  m/s),  respectively.  These 

figures  show  that  approximately  X sec  are  required  to  reach  the  steady-state  configuration. 


x(fT)  ,(FTI 

0 123456  0 12346678 


Figure  8 - Configuraiion  of  Nodes  at  Figure  9 - Configuration  of  Nodes  at 

Various  Times  during  Deployment.  Various  Times  during  Deployment. 

C = 1 Knot  C = 2 Knots 


I 

I 


8 


10 


11 


12 


f 


I 


Figure  10  — Configuration  of  Nodes  at  Various  Times 
during  Deployment.  C = 4 Knots 

In  .ill  ilircc  ca\cs.  the  final  deployed  configuration  modeled  by  the  four  nodes  v>as  in  close 
ipuciiu-nt  with  the  exact  configuration  calculated  by  using  the  differentia!  equations  (1 ) to 
i.st  li  IS  of  interest  to  note  that  because  of  its  large  inertia,  the  lower  weight  characteristically 
l.igs  during  the  initial  instants  of  deployment. 

l.ilde  2 shows  that  the  computer  execution  time  for  a total  of  20  sec  of  deployment 
tinie.  which  is  well  in  excess  of  the  8 sec  required  to  reach  the  steady-state  configuration,  is 
.ipprosiniately  100  sec. 


PROBLEM  4 - COMPLETE  MOORED  BUOY-CABLE-BODY  SYSTEM 

Piutilenr  Compute  the  dynamic  motions  for  the  moored  buoy-cable-body  system,  with 
[ui.iiiiclers  as  shown  in  Table 
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TABLE  7 - PARAMETERS  FOR  MOORED  BUOY-CABLE-BODY  SYSTEM 


4.8  ft  {1.46  nr>» 


Wind  Load  4.5  lb  (20.03  N) 


YCG  =3.7  ft  (1.13  m) 


Fluid  density  = 1.99  slugs/ft^  (1026  kg/m^) 


Lg  = 900  ft  (274.5  m) 


Cable  1 


Body  1 


Spheroidal  Surface  Buoy 
Waterplane  radius 
Draft 

Weight  in  air 
Moment  of  inertia 
Submerged  x-drag  area 


0.53  ft  (0.162  m) 

5.3  ft  (1.62  m) 

38  lb  (169.1  N) 

1.8  slug  ft^  (2.44  kg  m^) 
3.15  ft^  (0.293  m^; 


1500  ft  (457.5  m) 


L„  = 900  ft  (274.5  m) y 


Cable  1 


Body  2 I 

L„  = 500ft 
(152.5  m)  / 


Submerged  y-drag  area 

0 22  ft^  (0.020 

m*) 

Body  1 

Body  2 

Body  3 

Weight  in  fluid,  lb  (4.45  N) 

-40.0 

10.0 

15.0 

x-drag  area,  ft^  (0.0929  m^) 

0.44 

0.18 

0.70 

y-drag  area,  ft^  (0.0929  m^) 

0.44 

1.95 

0.08 

x-virtual  mass,  slugs  (14.6  kg) 

0.73 

0.55 

1.92 

y-virtual  mass,  slugs  (14.6  kg) 

0.73 

2.75 

1.35 

Cable  2 


Body  3 


L„  = 200  ft 
(61.0  m)^ 


Cable  2 


Cable  1 Cable  2 
Diameter,  in.  (2.54  cm)  0.15  0.25 

Normal  drag  coefficient  1.4  1 .4 

Tangential  drag  coefficient  0.02  0.02 

Weight  in  fluid,  Ib/ft  ( 1 4.6  N/m)  0.04  0.02 

Mass,  slugs/ft  (47.84  kg/m)  0.0015  0.0013 

Reference  tension,  lb  (4.45  N)  0.  1(X). 

C,  lb  (4.45  N)  22.000.  1,600. 

C,  1.0  1.5 

Internal  damping  coefficient,  0.  0. 

lb  sec  (4.45  N sec) 


Current  Profile; 


Depth 

(ft) 

(0.305  m) 
0. 
500. 
1000. 
1500. 


Current 
(knots) 
(0.515  m/s) 
2.50 
1.30 
0.50 
0.50 
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The  upper  body  is  a buoyancy  bag,  the  middle  body  is  a damping  disk,  and  the  lower  body 
is  an  acoustic  unit.  Consider  the  following  three  cases; 

Case  A.  Use  the  formulation  contained  in  the  program  for  a spheroidal  buoy.  Let  the  surfao. 
wave  be  composed  of  a single  component  with  frequency  = 0.1  cps.  amplitude  = 7.5  ft 
(2.29  m),  and  phase  = 0 degrees. 

Case  B.  Same  as  Case  A except  that  the  buoy  is  read  in  as  a spar  buoy,  modeled  by  1 1 

cross-sectional  areas  as  a function  of  depth  from  0 to  5.3  ft  (1.62  m). 

Case  C.  Same  as  Case  A except  that  the  surface  wave  amplitudes  are  to  be  calculated  by 

using  the  Pierson-Mosk'^witz  sea  spectrum  for  six  components.  Take  the  significant  wave 
height  as  15  ft  (4.58  m),  the  range  of  frequencies  from  0.04  to  0.28  cps.  and  the  phase  of 
the  lowest  frequency  component  equal  to  -60  deg. 

For  all  three  cases,  model  the  cable  by  4 nodes,  with  three  of  them  corresponding  to 
the  three  intermediate  bodies  and  the  extra  node  200  ft  of  cable  below  the  surface  buoy. 
Compute  the  dynamic  motions  for  50  sec.  For  the  initial  10  sec.  print  out  the  dynamic 
motions  everj'  0.1  sec.  For  the  final  40  sec,  increase  the  printout  interval  to  0.5  sec.  Let 
the  system  start  from  rest  with  the  initial  angle  and  tension  of  each  segment  equal  to  the 
steady-state  values.  Let  the  initial  x-displacemcnt  of  the  buoy  center  of  gravity  be  7.5  ft 
(2.29  m>. 

Solution:  The  data  cards  for  this  problem  arc  listed  in  Table  8.  ,\gain,  cards  which  are  the 
same  for  all  the  cases  arc  listed  at  the  top  of  the  table.  Tlie  cards  for  the  title  surface  buoy, 
and  surface  waves,  which  differ  for  each  case,  are  listed  at  the  bottom  of  the  table. 

Figure  1 1 shows  the  pitch  angle  tg  of  the  surface  buoy  for  all  three  c.oes  for 
20  < t < 40.  Perhaps  the  principal  feature  is  the  companson  of  the  rcMih  lor  the  same 
surtacc  buoy  treated  as  a spheroidal  buoy  and  as  a spar  buoy.  The  figure  show-s  that  when 
the  spheroidal  buoy  formulation  was  used,  pitch  results  vere  about  1 to  2 deg  higher  than 
those  obtained  from  the  spar  buoy  formulation.  This  disca'p^'n'-y  i''  due  to  the  difference  in 
the  added  mertia  terms  for  the  two  formulations  (compare  Equations  (27l,  (2Xi.  (32)  and 
Table  1 ». 

The  figure  also  shows  that  results  for  the  single  frequency  cases  4.\  and  4’t  exhibit  a 
periodic  behavior  with  a period  of  10  sec;  the  results  for  the  nuiltifreunency  c.ise  4C  show  a 
iiiore  random  behavior. 

Table  2 shows  that  the  execution  times  for  the  singic-freijuency  cases  A.\  ,md  4B  were 
approximately  80  sec.  This  lime  was  increased  to  iOO  sec  for  Case  4C.  wiierc  the  program 
must  calculate  six  components  to  obtain  the  surface  wave. 


TABLE  8 - INPUT  DATA  FOR  SAMPLE  PROBLEM  4 


1 

11 

21 

31 

41 

51 

61 

Card  1 

bbl 

Card  11 

1.99 

0. 

4.5 

-100. 

3.15 

1.0 

1.0 

Card  12 

10. 

0.1 

50. 

0.5 

1.0 

0. 

99999. 

Card  13 

200. 

700. 

900. 

500. 

200. 

Card  14 

0.15 

0.15 

0.15 

0.25 

0.25 

Card  15 

1.4 

1.4 

1.4 

1.4 

1.4 

Card  16 

0.02 

0.02 

0.02 

0.02 

0.02 

Card  17 

0.04 

0.04 

0.04 

0.02 

0.02 

Card  18 

0.0015 

0.0015 

0.0015 

0.0013 

0.0013 

Card  19 

0. 

0. 

0. 

100. 

100. 

Card  20 

22000. 

22000. 

22000. 

1600. 

1600. 

Card  21 

1.0 

1.0 

1.0 

1.5 

1.5 

Card  22 

0. 

0. 

0. 

0. 

Card  23 

0. 

1 

p 

10. 

15. 

Card  24 

0. 

0.44 

0.18 

0.70 

Card  25 

0. 

0.44 

1.95 

0.08 

0. 

Card  26 

0. 

0.73 

0.55 

1.92 

0. 

Card  27 

0. 

0.73 

2.75 

1.35 

0. 

Card  28 

0. 

500. 

1000. 

1500. 

Card  29 

2.50 

1.30 

0.50 

0.50 

Card  30 

99^. 

0. 

0. 

0. 

0. 

Card  31 

0. 

0. 

0. 

0. 

0. 

Card  32 

0. 

0. 

0. 

0. 

0. 

Card  33 

0. 

0. 

0. 

0. 

0. 

Card  34 

0.22 

38. 

GO 

1 

0. 

1.6 

3.7 

1.8 

Card  35 

7.5 

0. 

999. 

0. 

0. 

0. 

Card  2 

bbb  PROBLEM  4A,  SPHEROIDAL  BUOY,  SINGLE  WAVE  FREOUENCY 

Card  3 

bb1bb1bb5bb4bb2 

Card  4 

3500. 

Card  5 

0.53 

Card  6 

5.3 

Card  7 

0. 

Card  8 

7.5 

Card  9 

0.1 

Card  10 

0. 

5‘) 


TABLE  8 - (continued) 


Card  2 
Card  3 
Card  4 
Card  4 


1 11  21  31  41  51 

bbb  PROBLEM  4B,  SPAR  BUOY,  SINGLE  WAVE  FREQUENCY 
b11bb1bb5bb4bb2 


Card  5 

0.883 

0.852 

0.812 

0.757 

1.687 

0.600 

0.498 

0.380 

Card  5 

0.247 

0.097 

0. 

Card  6 

0. 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

Card  6 

4.5 

5.0 

5.3 

Card  7 
Card  7 
Card  8 
Card  9 
Card  10 
Card  2 
Card  3 
Card  4 
Card  5 0.53 

Card  6 5.3 

Card  7 0. 

Card  8 1015 

Card  9 0.04 

Card  10  -60. 


bbb  PROBLEM  4C,  SPHEROIDAL  BUQV  SIX  WAVE  FREQUENCfCS 

bb1bb6bb5bb4bb2 

3500. 

0.53 


60 
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APPENDIX 

LISTING  OF  COMPUTER  PROGRAM 
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